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Abstract 

Ultra  wideband  (UWB)  signals  characteristically  occupy  a  very  large  band¬ 
width  resulting  from  extremely  short  duration  pulses.  Given  an  impulse-like  UWB 
signal  occurs  so  quickly  and  its  energy  is  spread  across  such  a  large  bandwidth,  it  is 
extremely  difficult  to  detect  without  having  prior  knowledge  of  its  existence. 

This  research  develops  and  evaluates  techniques  for  the  non-cooperative  (non- 
matched  filter)  detection  of  impulse-like  UWB  signals  using  channelized  receiver 
architectures.  Each  technique  considered  is  modeled  and  simulations  conducted  to 
characterize  detection  performance,  the  results  of  which  are  compared  with  the  de¬ 
tection  performance  of  three  receivers:  the  matched  filter  receiver,  which  provides 
optimum  detection  performance  in  AWGN;  the  radiometer,  or  energy  detector;  and 
the  multi-aperture  cross  correlation  receiver.  The  input  signal  considered  is  a  single 
UWB  impulse  in  AWGN. 

It  is  shown  that  a  channelized  receiver  (with  no  downconversion)  can  provide 
approximately  2.5  dB  improvement  over  the  radiometer  when  performing  detec¬ 
tion  using  the  temporal-temporal  matrix  (TTM).  The  TTM  processing  technique 
provides  the  best  performance  of  all  the  proposed  channelized  receiver  techniques. 
Although  cross  temporal  matrix  (CTM)  detection  performance  exceeds  that  of  the 
radiometer  as  well,  it  does  not  perform  quite  as  well  as  TTM  detection  for  all  chan¬ 
nel  bandwidths.  Two  spectral-based  techniques  were  considered  as  well  with  neither 
performing  as  well  as  the  radiometer. 

Detection  with  a  downconverting  channelized  receiver  is  shown  dependent  on 
mixer  phase  value  with  performance  variation  generally  minimized  as  the  number 
of  channels  increases  (channel  bandwidth  decreases).  Although  not  always  the  case, 
there  are  combinations  of  signal-to-noise  ratio  and  channel  bandwidth  whereby  the 
downconverting  channelized  receiver  performs  better  than  the  radiometer. 


NON-COOPERATIVE  DETECTION  OF  ULTRA  WIDEBAND 


SIGNALS 


1.  Introduction 

1.1  Introduction 

The  FCC’s  First  Report  and  Order  released  in  April  2002  made  available  the 
unlicensed  spectrum  between  3.1  and  10.6  GHz  for  ultra  wideband  (UWB)  sys¬ 
tems  [1],  an  action  which  generated  renewed  interest  in  impulse  signaling  systems 
and  research.  Various  people  and  groups  have  shown  that  operating  an  ultra  wide¬ 
band  system  in  the  vicinity  of  existing  systems,  including  radars  [2],  narrowband 
military  receivers  [3],  802.11  wireless  LANs  [4-6],  and  GSM900,  UMTS/WCDMA, 
and  GPS  [7],  can  have  a  detrimental  effect  on  the  performance  of  such  systems.  It 
would  be  beneficial  to  these  systems  to  have  the  ability  to  detect  the  presence  of  a 
UWB  signal  and  adjust  operation  characteristics  to  improve  performance,  if  possible. 

1.2  Problem  Statement 

Ultra  wideband  (UWB)  signals  characteristically  occupy  a  very  large  band¬ 
width  resulting  from  extremely  short  duration  pulses.  Given  an  impulse- like  UWB 
signal  occurs  so  quickly  and  its  energy  is  spread  across  such  a  large  bandwidth,  it  is 
extremely  difficult  to  detect  without  having  prior  knowledge  of  its  existence.  This  re¬ 
search  focuses  on  developing  techniques  for  the  non-cooperative  (non-matched  filter) 
detection  of  impulse-like  UWB  signals  using  channelized  receiver  architectures. 

1.3  Research  Assumptions 

The  following  assumptions  were  used  throughout  the  research: 
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•  The  channel  is  modeled  as  being  additive  white  Gaussian  noise  (AWGN). 

•  The  impulse-like  UWB  signal  is  the  only  signal  present  in  the  RF  environment. 

•  All  RF  filter  spectral  responses  are  centered  at  the  UWB  signal  center  fre¬ 
quency.  The  RF  filter  bandwidth  is  fixed  when  comparing  results  across  re¬ 
ceiver  types  and  any  signal-to-noise  ratio  (SNR)  measurements  are  made  at 
the  output  of  the  RF  filter. 

•  All  signals  are  received  along  a  line-of-sight  path  from  the  transmitter.  No 
multi-path  signals  are  present. 

•  All  signal  detection  is  performed  using  test  statistics  generated  under  constant 
false  alarm  rate  (CFAR)  conditions,  i.e. ,  the  detection  threshold  varies  as  a 
function  of  signal-to-noise  ratio  (SNR)  such  that  a  constant  probability  of  false 
alarm  of  Pfa  —  10~2  is  maintained. 

1.4  Research  Scope 

Channelized  receiver  processing  techniques  were  developed  for  the  purpose  of 
detecting  impulse-like  UWB  signals.  Each  technique  considered  is  modeled  and  sim¬ 
ulations  conducted  to  characterize  detection  performance,  the  results  of  which  are 
compared  with  the  detection  performance  of  three  receivers:  the  matched  filter  re¬ 
ceiver,  which  provides  optimum  detection  performance  (highest  probability  of  detec¬ 
tion  for  a  given  signal-to-noise  ratio)  in  AWGN;  the  radiometer,  or  energy  detector; 
and  the  multi-aperture  cross  correlation  receiver. 

1.5  Research  Approach 

Initial  research  into  impulse-like  UWB  signals  provides  analytic  derivation  of 
power  spectral  density  expressions  for  signals  employing  various  UWB  modulations, 
including  pulse  position  modulation  (PPM),  pulse  amplitude  modulation  (PAM), 
biorthogonal  pulse  position  modulation  (BPPM),  and  a  uniform  pulse  train  (UPT) 
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as  commonly  used  in  radar  applications  (no  random  data  modulation).  The  analytic 
results  for  each  modulation  type  are  then  plotted  to  illustrate  the  characteristics  of 
each  spectrum. 

A  matched  filter  receiver  model  is  developed  next  based  on  the  correlator  imple¬ 
mentation  of  the  matched  filter.  This  receiver  model  assumes  the  shape  and  arrival 
time  of  the  received  UWB  pulse  are  known  a  priori.  Thus,  the  matched  filter  receiver 
provides  optimum  detection  performance  in  AWGN  and  results  for  this  model  serve 
as  an  upper  bound  (best  case)  on  achievable  detection  performance  for  the  proposed 
channelized  detection  methods.  Two  non-cooperative  receivers  are  also  considered, 
including  the  radiometer  and  multi-aperture  cross  correlation  receiver.  Simulated 
detection  results  for  these  receivers  serve  as  a  baseline  for  detection  performance 
using  proposed  channelized  methods. 

The  research  concludes  with  the  development  and  detection  performance  char¬ 
acterization  of  a  channelized  receiver  and  downconverting  channelized  receiver  model. 
Both  the  channelized  and  downconverting  channelized  receiver  spectrally  divide 
(channelize)  the  input  signal  into  channels  using  bandpass  filters.  However,  the 
downconverting  channelized  receiver  also  employs  mixers  and  lowpass  filters  to  inde¬ 
pendently  downconvert  each  channelized  signal  to  baseband.  For  both  receivers,  the 
final  filtered  signals  are  digitized  in  preparation  for  digital  processing.  The  digital 
processor  uses  data  from  a  finite  length  observation  interval  which  is  stored  in  ma¬ 
trices  of  dimension  M  x  S,  where  M  is  the  number  of  receiver  channels  and  S  is  the 
number  of  time  samples  spanning  the  interval.  The  detection  performance  of  each 
processing  technique  is  determined  through  simulation  and  subsequently  compared 
with  radiometric  and  matched  filter  results. 

1.6  Materials  and  Equipment 

The  receiver  architectures  and  detection  techniques  presented  in  this  work 
were  simulated  using  MATLAB®  Version  7.0,  developed  by  Mathworks,  Inc.  The 
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simulations  were  run  on  a  2.2  GHz  Athlon  XP  personal  computer  and  an  AFIT 
mainframe  UNIX  based  system. 

1.7  Thesis  Organization 

Chapter  2  provides  background  information  on  the  impulse-like  ultra  wideband 
signals,  including  different  modulation  types  and  their  associated  power  spectral 
densities.  The  receiver  architectures  used  for  both  cooperative  and  non-cooperative 
signal  detection  are  introduced  as  well.  Chapter  3  presents  the  methodology  used 
for  conducting  the  research,  including  a  detailed  discussion  of  how  the  four  different 
data  matrices  (two  time-based  and  two  spectral-based)  are  generated  from  channel¬ 
ized  receiver  processing.  The  chapter  concludes  with  a  description  of  channelized 
receiver  detection  using  each  of  these  four  matrices.  Chapter  4  provides  simulated 
detection  results  and  analysis  for  both  the  cooperative  and  non-cooperative  receivers 
considered.  Chapter  5  presents  conclusions  drawn  from  the  research  and  provides  rec¬ 
ommendations  for  possible  future  research.  Additional  supporting  data  is  provided  in 
the  appendices:  Appendix  A  contains  simulated  detection  results  for  receiver  chan¬ 
nel  bandwidths  not  included  in  Chapter  4;  Appendix  B  contains  processed  matrix 
data  and  simulation  results  for  an  impulse-like  UWB  signal  occurring  at  a  different 
arrival  time  than  what  was  used  for  results  presented  in  Chapter  3  and  Chapter  4; 
Appendix  C  contains  the  MATLAB®  code  used  in  the  simulations. 
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2.  Background 


2.1  Introduction 

This  chapter  introduces  the  characteristics  of  various  ultra  wideband  (UWB) 
signals  and  the  receivers  that  will  be  used  for  UWB  signal  detection.  Section  2.2  de¬ 
scribes  various  UWB  modulations  and  associated  power  spectral  densities.  Section 
2.3  describes  matched  filter  detection,  and  Section  2.4  describes  the  non-cooperative 
detection  receivers  considered  for  this  research,  including  the  radiometer,  the  multi- 
aperture  cross  correlation  receiver,  and  the  channelized  receiver.  Section  2.5  sum¬ 
marizes  the  chapter. 

2.2  Ultra  Wideband  Overview 

Ultra  wideband  (UWB)  systems  transmit  data  over  a  very  large  bandwidth 
using  extremely  short  duration  pulses.  The  FGC’s  First  Report  and  Order  released 
in  April  2002  [1]  places  emissions  regulations  on  the  effective  isotropic  radiated  power 
(EIRP)  of  UWB  devices  to  limit  potential  interference  with  narrowband  systems 
utilizing  the  same  spectrum.  Specifically,  the  FCC  specifies  a  maximum  allowable 
UWB  EIRP  of  41.3  dBm  in  the  3.1  to  10.6  GHz  band  (measured  at  -10  dB  points) 
for  indoor  and  unrestricted  handheld  devices. 

2.2.1  Ultra  Wideband  Definition.  According  to  the  FCC  [1],  a  system 
is  considered  UWB  if  it  has  1)  a  bandwidth  greater  than  or  equal  to  500  MHz,  or 
2)  a  fractional  bandwidth  greater  than  20%.  Fractional  bandwidth  Bj  is  given  by 

B’  =  2{% rf)  (2J) 

where  /l  and  fH  are  the  lower  and  upper  —10  dB  bandwidth  frequencies  of  a  signal. 
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2.2.2  Ultra  Wideband  Communication  Signaling.  The  analytic 
form  that  is  often  chosen  to  represent  the  transmitted  UWB  pulse  is  the  Gaussian 
monocycle.  Accounting  for  both  transmit  and  receive  antenna  effects,  the  received 
UWB  waveform  is  the  second  derivative  of  a  Gaussian  impulse  and  is  given  by 


w(t) 


(  t  \2 

\  ftV  1 

1  -  4tt  — 

exp 

-2tt  — 

V  Tm  )  _ 

\TmJ 

(2.2) 


where  the  impulse  width  parameter  rm  is  approximately  0.4  times  the  pulse  width 
Tw  [8].  A  single  UWB  pulse  is  plotted  in  Fig.  2.1  using  Tw  —  0.4  ns. 


Figure  2.1  Received  Gaussian  Impulse  Accounting  for  Transmit /Receive  Antenna 
Effects 

The  general  expression  for  a  received  UWB  signal  using  an  arbitrary  modula¬ 
tion  is 


OO 

s(t)  =  ^2  Ak-w(t-  kTs  -  Bk A)  (2.3) 

k= — oo 
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where  Ak  and  Bk  are  specified  according  to  modulation  type,  Ts  is  the  symbol  dura¬ 
tion,  and  A  is  the  relative  position  modulation  offset.  Four  types  of  UWB  signals  are 
considered  in  this  work  by  specifying  various  values  for  Ak  and  Bk  in  (2.3),  including: 

•  Uniform  Pulse  Train  (UPT):  Ak  —  Bk  —  1 

•  Pulse  Position  Modulation  (PPM):  Ak  —  1  and  Bk  —  ±1 

•  Pulse  Amplitude  Modulation  (PAM):  Ak  —  ±1  and  Bk  —  1 

•  Biorthogonal  Pulse  Position  Modulation  (BPPM):  Ak  —  ±1  and  Bk  —  ±1 

As  shown,  the  PPM  and  PAM  signals  are  binary  data  modulated  waveforms,  i.e., 
either  Ak  or  Bk  has  the  ±  sign  variation  for  a  given  modulation  type  with  one  bit 
value  determining  the  appropriate  amplitude  sign.  The  BPPM  signal  is  a  4-ary  data 
modulated  waveform  and  thus  both  Ak  and  Bk  terms  have  the  ±  sign  variation  with 
two  bit  values  (one  per  term)  determining  the  amplitude  signs. 

2.2.3  Ultra  Wideband  Power  Spectral  Density.  The  general  analytic 
expression  for  the  power  spectral  density  (PSD)  of  an  ultra  wideband  signal  is  derived 
using  a  method  similar  to  those  presented  in  [9-11].  To  simplify  the  analysis,  the 
signal  at  the  receiver  given  in  (2.3)  will  be  described  in  terms  of  the  modulation 
process  and  the  pulse  shape  as  s(t)  —  m(t )  *  w(t),  where  *  is  the  convolution 
operation  and  m(t),  the  modulation  process,  is  given  by: 

OO 

m(t)  =  ^  Ak  ■  5(t  —  kTs  -  Bk A)  (2.4) 

k=— oo 

The  modulation  process  PSD  is  found  by  first  calculating  its  autocorrelation 
function,  i?mof|(r).  It  must  be  noted  that  Rmod{ t)  is  periodic  in  r  with  Ts  and  the 
time-average  autocorrelation  function  must  be  derived  [10]. 

Rmod(r)  =  [  E  [m  (t)  m(t  +  r)]  dt  (2.5) 

hi  0 
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^orf(r)  =  y^E  [A°Al]  E\-6^t-  lT°  ~  ~  B °)  A)]  (2'7) 

s  l=—oo 

The  above  expressions  assume  Ak  and  Bk  are  uncorrelated.  The  time-average  auto¬ 
correlation  is  next  split  into  two  parts:  l  —  0  and  l  ^  0. 


Rmodir)  =  ^l<J(r)  +  E  E  [^]  ■  E  [6  (t-  ITS  -  (Bt  -  B0)  A)]  (2.8) 

s  s  i^o 

The  modulation  process  PSD,  Smod(f),  is  then  found  by  taking  the  Fourier  transform 
of  the  time-average  autocorrelation  function  in  (2.8)  . 


Smod(f)  =  £  ElAoAfrEle-R^-^]- 

-L  s  -L  s  , 

/=— oo 

Once  the  modulation  process  PSD  has  been  found,  the  total  UWB  signal  PSD  is 
given  by: 


-<*(/-  — 
T„  V  T, 


-  1 


(2.9) 


SuwB(f)  =  Smod(f)\W(f)\- 


(2.10) 


where  W(f )  is  the  Fourier  transform  of  the  received  UWB  waveform,  w(t),  given  in 
(2.2).  It  can  be  shown  that  the  Fourier  transform  of  w(t)  is 


7 TT^ 

W{f)  =  •  exp 
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2.2.3. 1  PSD  of  Uniform  Pulse  Train.  A  uniform  train  of  UWB 
pulses  is  created  using  Ak  —  1  and  Dk  —  1  in  (2.3).  Using  these  values,  the  ana¬ 
lytic  expression  for  the  UPT  power  spectral  density  can  be  found  by  evaluating  the 
expectations  in  (2.9)  which  are  given  by: 

E  [Aq]  =  E  [A0Ai]  =  1 

E  reJ'27r/(B!-Bo)Aj  =  1 


The  PSD  for  the  modulation  process  then  has  the  following  analytic  expression: 


1  /  / 
Suprtf)  =  7fv  6  (  f 


l=— OO 


Tx 


(2.12) 


and  the  total  UPT  power  spectral  density  is  found  by  substituting  (2.12)  into  (2.10) 
for  Smod(f). 


2. 2. 3. 2  PSD  of  Binary  PPM  Signal.  A  binary  pulse  position 
modulated  (PPM)  UWB  signal  is  created  using  Ak  —  1  and  allowing  Bk  =  ±1  in 
(2.3).  Using  P[Bk  —  1]  =  P[Bk  —  —  1]  =  0.5  and  evaluating  the  expectations  in  (2.9) 
gives: 

E  [Aq]  =  E  [A0A;]  =  1 
E  [ej27r^^Bi_Bo^A]  =  -  [1  +  cos  (47t/A)] 

This  yields  a  PSD  for  the  modulation  process  given  by 

SppM(f)  =  2jT  [  1  -  cos  (4tt/A)]  +  ^  [1  +  cos  (4tt/A)] 

The  total  PSD  is  found  by  substituting  (2.13)  into  (2.10)  for  Smod{f ) 


4)  <213) 
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2. 2. 3. 3  PSD  of  Binary  PAM  Signal.  A  binary  pulse  amplitude 
modulated  (PAM)  UWB  signal  is  created  by  allowing  —  ±1  and  setting  —  1 
in  (2.3).  Using  P[Ak  —  1]  =  P[Ak  —  —1]  =  0.5,  the  expectations  in  (2.9)  are  given 
by: 


E  K]  =  1 

E[A0At]  =  0,  VZ^O 

E  rej27T/(Bi-s0)Aj  =  1 

This  yields  a  PSD  for  the  modulation  process  given  by 

Spam(/)  =  ±  (2-14) 

-L  s 

The  total  PSD  is  found  by  substituting  (2.14)  into  (2.10)  for  Smod{f)- 

2. 2. 3. 4  PSD  of  4-Ary  Biorthogonal  PPM  Signal.  A  4-ary 
biorthogonal  pulse  position  modulated  (BPPM)  signal  combines  the  characteristics 
of  the  binary  PPM  and  PAM  signals  to  create  a  4-ary  signaling  scheme.  A  BPPM 
signal  is  created  by  allowing  A^  —  ±1  and  Bk  —  ±1  in  (2.3).  Using  P[Ak  =  1]  = 
P[Ak  =  —  1]  =  0.5  and  P[Bk  =  1]  =  P[Bk  =  —  1]  =  0.5,  the  analytic  expression  for 
the  PSD  can  be  found  after  evaluating  the  expectations  in  (2.9). 

E  [Ag]  =  1 
E[A<)Al]  =  0,  V/^0 
E  =  1  [i  +  CQS  (47r/A)] 

This  yields  a  PSD  for  the  modulation  process  given  by 

SbppmU )  =  Tfp  (2.15) 

■L  s 
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The  total  PSD  is  found  by  substituting  (2.15)  into  (2.10)  for  Smod(f ). 

2. 2. 3. 5  Comparison  of  UWB  Signal  PSDs.  Analytic  results 
derived  above  are  shown  plotted  in  Fig.  2.2  using  parameter  values  of  Tw  —  0.4  ns, 
A  =  Tw/2  —  0.2  ns,  and  rm  =  0.4 Tw  —  0.16  ns.  Both  PPM  and  BPPM  use  a  symbol 
duration  of  Ts  —  2 Tw  and  PAM  uses  a  symbol  duration  of  Ts  —  Tw.  The  uniform 
pulse  train  uses  a  pulse  repetition  interval  (PRI)  of  Ts  —  2 Tw. 


Figure  2.2  Analytic  PSDs  of  Various  UWB  Waveforms 


2.3  Cooperative  Detection 

2.3.1  Matched  Filter  Detection.  A  receiver  can  achieve  optimum 
detection  performance  in  additive  white  gaussian  noise  (AWGN)  if  the  form  and 
arrival  time  of  the  received  UWB  pulse  are  known.  When  these  parameters  are 
known,  the  correlator  implementation  of  a  matched  filter  (MF)  receiver  as  shown  in 
Fig.  2.3  [12]  can  be  used  to  detect  and  estimate  signal-of-interest  (SOI)  presence.  In 
Fig.  2.3,  r(t )  is  the  received  signal  comprised  of  a  desired  SOI  component  s(t)  and 
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an  AWGN  component  n(t).  Although  not  functionally  part  of  the  MF  process,  RF 
filter  Wmf  has  been  incorporated  for  consistency  with  subsequent  detectors  to  be 
considered.  The  detection  process  is  completed  by  comparing  the  MF  output  test 
statistic  Z  with  a  threshold.  If  Z  exceeds  the  threshold  signal  presence  is  declared. 

s(t) 


Figure  2.3  Matched  Filter  Detection  Implemented  with  Correlation 


2.4  Non-Cooperative  Detection 

Assuming  signals  in  the  environment  are  unknown,  non-matched  filter  detec¬ 
tion  techniques  must  be  employed.  Non-matched  filter  detection  is  perhaps  the  most 
challenging  detection  approach  that  can  be  taken,  and  the  problem  is  even  more 
challenging  when  multiple  signals  exist  simultaneously. 

As  with  the  matched  filter  detector  described  in  the  previous  section,  all  non- 
cooperative  detection  techniques  considered  in  this  work  generate  a  test  statistic  Z 
for  threshold  comparison.  The  test  statistic  output  from  each  processing  technique  is 
compared  with  threshold  ZT  which  can  be  established  using  various  detection  criteria 
(e.g.  Bayes,  Neyman-Pearson,  etc.)  [13].  The  SOI  is  declared  “present”  whenever  Z 
exceeds  ZT.  Assuming  the  SOI  is  present  during  test  statistic  generation,  one  of  two 
conditions  occurs,  either  1)  Z  exceeds  the  threshold,  valid  detection  occurs,  and  a 
probability  of  detection  (Pd)  can  be  computed,  or  2)  Z  falls  below  the  threshold,  no 
detection  occurs,  and  a  probability  of  miss  (Pm)  can  be  computed.  When  the  SOI 
is  not  present  during  processing,  only  channel  noise  is  present  and  two  additional 
conditions  are  possible,  either  1)  Z  exceeds  the  threshold,  in  which  case  a  false 
detection  occurs  and  a  probability  of  false  alarm  (Pda)  can  be  computed,  or  2)  Z  falls 
below  the  threshold  and  no  detection  occurs  as  desired.  The  threshold  can  adapt  to 
changing  channel  conditions  such  that  a  constant  Pfa  is  maintained.  This  is  known 
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as  constant  false  alarm  rate  (CFAR)  processing  [14].  The  following  sections  introduce 
three  receiver  structures  that  are  considered  for  non-cooperative  detection  as  part 
of  this  work,  including  the  radiometer,  multi-aperture  cross  correlation  receiver,  and 
channelized  receiver. 


2.4.1  Radiometric  Detection.  As  shown  in  the  radiometer  block  dia¬ 
gram  of  Fig.  2.4,  a  radiometer  detects  signal  energy  in  bandwidth  Wra d  Hz  using  a 
coherent  processing  time  of  T  sec  [15].  The  resultant  test  statistic  Z  is  compared  to 
threshold  value  ZT,  if  Z  >  ZT  signal  presence  is  declared  (detection)  and  if  Z  <  ZT 
detection  does  not  occur. 


>  =>  Signal  Present 
<  =4>  Not  Present 


Figure  2.4  Radiometric  Detection  Using  Bandwidth  W rad  and  Coherent  Integra¬ 
tion  Time  T 


2.4.2  Multi- Aperture  Cross  Correlation  Receiver  Detection.  The 

multi-aperture  cross  correlation  receiver  in  Fig.  2.5  intercepts  a  signal  at  two  spa¬ 
tially  separated  antennas,  filters  the  received  signals,  calculates  the  cross  correlation 
(denoted  by  ®)  for  various  relative  arrival  times,  and  compares  the  peak  output  re¬ 
sponse  to  a  threshold  set  using  Bayes,  Neyman-Pearson,  or  other  detection  criteria. 
If  the  peak  value  exceeds  the  threshold,  a  signal  is  declared  present.  This  implemen¬ 
tation  is  an  approximation  of  the  correlator  implementation  of  the  matched  filter 
receiver  in  Fig.  2.3  [16].  The  major  difference  here  is  that  instead  of  correlating  a 
noisy  received  signal  with  a  locally  generated  noise-free  signal  (as  done  with  a  refer¬ 
ence  signal  in  the  matched  filter  receiver),  two  independently  received  noisy  signals 
are  cross  correlated.  In  this  case,  the  relative  arrival  time  of  the  signals  need  not 
be  calculated  and/or  estimated  and  time  alignment  (synchronization)  with  a  locally 
generated  reference  is  not  required. 
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Figure  2.5  Multi-Aperture  Cross-Correlation  Receiver  Detection 

2.4.3  Channelized  Receiver  Detection.  Channelized  receiver  tech¬ 
niques  provide  perhaps  the  best  alternative  to  wideband  radiometric  detection,  pro¬ 
viding  wide  instantaneous  bandwidth  such  that  ah  signals  present  are  simultaneously 
received  without  tuning  the  receiver  to  a  specific  signal  [17].  Channelized  receivers 
can  be  implemented  using  either  analog  or  digital  techniques.  Regardless  of  the 
implementation,  the  fundamental  channelized  receiver  consists  of  a  bank  of  filters 
as  shown  in  Fig.  2.6,  the  outputs  of  which  are  processed  to  arrive  at  some  desired 
conclusion. 


Figure  2.6  UWB  Channelized  Receiver 


Unfortunately,  implementing  this  receiver  in  digital  hardware  may  not  be  prac¬ 
tical  since  the  UWB  frequency  range  extends  to  frequencies  higher  than  what  current 
analog-to-digital  converter  technology  supports  (i.e.  sampling  at  or  above  Nyquist 
is  not  feasible).  One  possible  solution  to  this  limitation  is  to  use  a  downconverting 
channelized  receiver  as  shown  in  Fig.  2.7,  an  approach  which  has  been  proposed 
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in  [18].  As  shown,  the  downconverting  receiver  uses  M  mixers  operating  at  equally 
spaced  frequencies  (/i,  f2,  .../m)  and  phase  values  (6\,  d2, ... 9m )  which  can  be  varied 
either  dependency  or  independently  as  a  function  of  implementation.  The  mixers 
are  followed  by  a  bank  of  M  lowpass  filters  to  decompose  the  received  signal  into  M 
subbands.  The  downconverted  outputs  are  then  processed  to  arrive  at  some  desired 
conclusion. 


/iA 


Figure  2.7  UWB  Downconverting  Channelized  Receiver 


2.5  Summary 

This  chapter  introduced  the  characteristics  of  an  ultra  wideband  signal,  in¬ 
cluding  the  various  modulation  types  and  associated  power  spectral  densities.  The 
receiver  structures  used  for  detection  in  this  work  were  also  introduced,  includ¬ 
ing  the  matched  filter,  radiometer,  multi-aperture  cross  correlation  receiver,  and 
channelized  receiver.  Operations  performed  by  the  matched  filter,  radiometer,  and 
multi-aperture  cross  correlation  receiver  are  known,  thus  the  detection  performance 
of  each  can  be  readily  analyzed  through  simulation.  Basic  channelized  receiver  struc¬ 
ture  and  operation  is  also  known.  However,  specific  digital  processing  techniques  for 
non-cooperative  detection  were  unavailable.  Therefore,  several  data  processing  tech¬ 
niques  are  introduced  in  Chapter  3  for  signal  detection. 
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3.  Methodology 


3.1  Introduction 

This  chapter  introduces  the  channelized  receiver  processing  methods  developed 
to  detect  UWB  signals.  Fundamentally,  the  goal  is  to  achieve  the  best  possible  de¬ 
tection  performance,  as  bounded  below  by  the  radiometer  and  above  by  the  matched 
filter.  Each  of  the  four  processing  techniques  considered  are  presented  in  Section  3.2. 
The  method  used  for  determining  probability  of  detection  under  fixed,  constant  false 
alarm  rate  conditions  is  presented  in  Section  3.3. 

3.2  Channelized  Receiver  Processing 

The  channelized  receiver  outputs  of  Figs.  2.6  and  2.7  are  used  to  form  an 
M  x  S  Channelized  Data  Matrix  (CDM),  where  M  is  the  number  of  channels  and 
S  is  the  number  of  time  samples  in  the  observation  interval.  This  Channelized  Data 
Matrix  is  processed  using  techniques  introduced  in  the  following  sections  with  focus 
on  achieving  the  best  possible  detection  performance. 

Parameters  for  generating  data  and  figures  presented  in  Chapters  3  and  4  are 
as  follows: 

•  UWB  Pulse  Width  -  0.4  ns 

•  Total  Signal  Duration  (Observation  Interval)  -  5.2  ns 

•  Time  Sample  Spacing  -  0.01  ns 

•  RF  Filter  Frequency  Range  (-4.0  dB  Bandwidth)  -  2.5  to  7.5  GHz 

•  Channel  Bandwidth  (-3.0  dB  Bandwidth)  -  250  MHz 

•  Number  of  Channels  (M)  -  20 

•  Downconversion  Mixer  Phase  Values  (at  start  of  observation  interval)  -  0° 
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A  plot  of  the  received  UWB  signal  based  on  these  parameters  is  shown  in  Fig.  3.1. 
Additional  data  and  figures  for  a  received  UWB  signal  time-offset  from  center  are 
given  in  Appendix  B. 


Figure  3.1  UWB  Channelized  Receiver  Input  Signal 

3.2.1  Temporal-Temporal  Matrix  (TTM).  The  temporal-temporal 
matrix  (TTM)  is  formed  by  performing  an  Appoint  inverse  Fourier  transform  (IFFT) 
on  each  column  (frequency  samples)  of  the  CDM  and  taking  the  absolute  value  of 
each  element  in  the  resultant  matrix  (to  eliminate  complex  terms) .  Zero-padding  is 
employed  if  there  are  more  IFFT  points  than  receiver  channels  (N  >  M)  and  trun¬ 
cation  is  employed  if  there  are  fewer  IFFT  points  than  receiver  channels  ( N  <  M). 
The  truncation  is  performed  such  that  samples  from  lower  frequency  channels  are 
maintained  and  samples  from  higher  frequency  channels  are  discarded  as  required. 
Representative  TTMs  for  the  received  UWB  pulse  of  Fig.  3.1  input  to  both  the  chan¬ 
nelized  receiver  and  the  downconverting  channelized  receiver  are  shown  in  Fig.  3.2 
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and  Fig.  3.3,  respectively.  These  matrices  were  formed  by  processing  the  CDM 
columns  using  64-point  IFFTs. 
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Figure  3.2  Temporal-Temporal  Matrix  (TTM)  for  Channelized  Receiver  Data,  No 
Noise  Present 

3.2.2  Cross  Temporal  Matrix  (CTM).  The  cross  temporal  matrix 
(CTM)  is  formed  by  calculating  the  absolute  value  of  the  correlation  between  all 
combinations  of  TTM  columns.  Individual  CTM  elements  are  given  by 

TTM 

where  Nt f  jt  is  the  number  of  IFFT  points  used  to  form  the  TTM.  The  entire  matrix 
can  be  constructed  more  efficiently  using  matrix  operations  per  the  following 

CTM  =  -J—  \TTMh  •  TTM  I  (3.2) 

Nifft 


('I'M(i.j)  Y, 


if  ft  rn_1 


*(m,  i)  ■  TTM(m,j)  ,  *  =  1, ...,  S,  j  =  1, ....  S  (3.1) 


Time 
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Figure  3.3  Temporal-Temporal  Matrix  (TTM)  for  Downconverting  Channelized 
Receiver  Data,  No  Noise  Present 

where  (■ )H  is  the  hermitian,  or  conjugate-transpose,  operation.  Representive  CTMs 
for  the  received  UWB  pulse  of  Fig.  3.1  input  to  both  the  channelized  receiver  and  the 
downconverting  channelized  receiver  are  shown  in  Fig.  3.4  and  Fig.  3.5,  respectively. 
These  matrices  were  formed  using  the  TTM  matrices  of  Fig.  3.2  and  Fig.  3.3. 

3.2.3  Spectral-Spectral  Matrix  (SSM).  The  spectral-spectral  matrix 
(SSM)  is  formed  by  performing  an  IV-point  Fourier  transform  (FFT)  on  each  CDM 
row  and  taking  the  absolute  value  of  each  element  in  the  matrix  (to  eliminate  com¬ 
plex  terms).  The  number  of  points  in  the  FFT  ( N ),  must  be  greater  than  or  equal 
to  the  number  of  samples  in  the  observation  interval  (S)  so  that  no  time  samples  are 
truncated  when  performing  the  FFT.  This  ensures  that  time  samples  corresponding 
to  a  received  UWB  pulse  are  not  inadvertantly  discarded  by  the  FFT  operation. 
Zero-padding  is  employed  if  there  are  more  FFT  points  than  samples  in  the  observa¬ 
tion  interval  (N  >  S ).  Representative  SSMs  for  the  received  UWB  pulse  of  Fig.  3.1 
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Figure  3.4  Cross  Temporal  Matrix  (CTM)  for  Channelized  Receiver  Data,  No 
Noise  Present 

input  to  both  the  channelized  receiver  and  the  downconverting  channelized  receiver 
are  shown  in  Fig.  3.6  and  Fig.  3.7,  respectively.  These  matrices  were  formed  by 
processing  the  CDM  using  512-point  FFTs. 


3.2.4  Cross  Spectral  Matrix  (CSM).  The  cross  spectral  matrix  (CSM) 
is  formed  by  calculating  the  absolute  value  of  the  correlation  between  all  combina¬ 
tions  of  SSM  columns.  Individual  CSM  elements  are  given  by 


CSM(i,j)  = 


M 


M 


Y  SSM*(m,  i)  ■  SSM(m,j) 


m—  1 


i  =  1  ,...,Nfft,  j  =  l,...,Nfft 


(3.3) 

where  Nfft  is  the  number  of  FFT  points  used  to  form  the  SSM.  The  entire  CSM  can 
be  constructed  more  efficiently  using  matrix  operations  per  the  following 


CSM  =  -1-  I SSMH  ■  SSM 
M  1 


(3.4) 
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Figure  3.5  Cross  Temporal  Matrix  (CTM)  for  Downconverting  Channelized  Re¬ 
ceiver  Data,  No  Noise  Present 

where  (• )H  is  the  hermitian,  or  conjugate-transpose,  operation.  Representative  CSMs 
for  the  received  UWB  pulse  of  Fig.  3.1  input  to  both  the  channelized  receiver  and  the 
downconverting  channelized  receiver  are  shown  in  Fig.  3.8  and  Fig.  3.9,  respectively. 
These  matrices  were  formed  using  the  SSM  matrices  of  Fig.  3.6  and  Fig.  3.7. 
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Figure  3.6  Spectral-Spectral  Matrix  (SSM)  for  Channelized  Receiver  Data,  No 
Noise  Present 


Frequency 


Figure  3.7  Spectral-Spectral  Matrix  for  Downconverting  Channelized  Receiver 
Data,  No  Noise  Present 
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Figure  3.8  Cross  Spectral  Matrix  (CSM)  for  Channelized  Receiver  Data,  No  Noise 
Present 
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Figure  3.9  Cross  Spectral  Matrix  (CSM)  for  Downconverting  Channelized  Re¬ 
ceiver  Data,  No  Noise  Present 


3.3  Channelized  Receiver  Detection 


Two-dimensional  threshold  detection  is  performed  on  each  of  the  processed 
matrices  (TTM,  CTM,  SSM,  CSM)  described  above  to  determine  probability  of  de¬ 
tection  (Pd)  for  that  particular  processing  method.  Note  that  for  both  the  TTM  and 
SSM  cases,  only  linear  transformation  (IFFT  and  FFT,  respectively)  of  the  channel¬ 
ized  data  matrix  occurs  and  there  is  no  anticipated  processing  gain  (improvement) 
in  terms  of  enhanced  detection  performance.  However,  both  the  CTM  and  CSM 
cases  involve  correlation  and  thus  some  processing  gain  is  anticipated. 

For  threshold  detection,  the  maximum  value  of  the  processed  matrix  is  chosen 
as  the  test  statistic  for  subsequent  comparison  with  the  threshold.  The  threshold 
used  in  the  detection  process  is  determined  by  processing  a  series  of  matrices  created 
using  R  noise  realization  inputs  (with  noise  power  set  according  to  desired  SNR)  and 
generating  a  set  of  R  test  statistics.  The  threshold  is  then  set  to  achieve  the  desired 
probability  of  false  alarm,  Pfa,  using  the  R  noise-only  test  statistics.  Note  that 
for  all  results  presented  in  this  work,  Pfa  is  maintained  constant  as  SNR  is  varied. 
This  condition  is  referred  to  as  constant  false  alarm  rate  (CFAR)  detection  [14]. 
Next,  matrix  data  is  regenerated  using  the  same  noise  realizations  with  the  signal 
of  interest  present.  Test  statistics  from  these  results  are  compared  to  the  threshold 
and  Pd  is  determined  as  the  number  of  times  the  threshold  is  exceeded  divided  by 
R  total  realizations.  This  process  is  repeated  for  each  desired  SNR  value. 

A  similar  detection  process  is  used  to  calculate  Pd  for  the  radiometer,  matched 
filter,  and  multi-aperture  cross  correlation  receiver  implementations,  although  the 
outputs  of  these  receivers  are  either  multi-valued  one-dimensional  signals  (multi¬ 
aperture  cross  correlation  receiver),  in  which  case  the  test  statistic  is  the  maximum 
value,  or  single  valued  (radiometer  and  matched  filter),  in  which  case  the  test  statistic 
is  the  output  itself. 
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Figure  3.10  illustrates  the  process  flow  used  for  channelized  receiver  detection. 
Various  parameters  are  changed  to  generate  results  for  different  detection  scenarios, 
including: 

•  Received  Signal-to-Noise  Ratio  (SNR)  (measured  at  the  output  of  the  RF  filter) 

•  Number  of  IFFT  and  FFT  points  used  to  form  the  TTM  (CTM)  and  SSM 
(CSM),  respectively 

•  Channelized  receiver  channel  bandwidths 

•  Starting  phases  of  downconversion  mixers  (downconverting  channelized  re¬ 
ceiver  only) 

3.4  Summary 

This  chapter  introduced  four  different  channelized  receiver  processing  tech¬ 
niques  considered  for  this  work,  including  the  temporal-temporal  matrix  (TTM), 
the  cross  temporal  matrix  (CTM),  the  spectral-spectral  matrix  (SSM),  and  the  cross 
spectral  matrix  (CSM).  A  method  for  determining  probability  of  detection  (Pd)  us¬ 
ing  threshold  detection  under  constant  false  alarm  rate  (CFAR)  conditions  was  also 
introduced.  Results  for  this  detection  process  with  each  of  the  matrices  generated 
are  shown  in  Chapter  4. 
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Channelized  Receiver  Filter  Output  Data  /  Matrix 
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Figure  3.10  Detection  Processing  Flow 
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4.  Detection  Results  and  Analysis 

4.1  Introduction 

This  chapter  provides  detection  performance  results  for  the  receivers  and  process¬ 
ing  techniques  introduced  in  Chapters  2  and  3.  First,  Section  4.2  provides  matched 
filter  and  radiometric  detection  performance  results  which  bound  the  desired  results 
for  the  proposed  techniques.  Next,  Section  4.3  presents  multi-aperture  cross  cor¬ 
relation  receiver  performance.  Finally,  Section  4.4  provides  detection  performance 
results  using  the  matrix  data  and  techniques  introduced  in  Chapter  3  for  both  the 
channelized  receiver  and  the  downconverting  channelized  receiver. 

4.2  Matched  Filter  and  Radiometric  Detection  Performance 

Matched  filter  and  radiometric  detection  results  are  used  as  relative  metrics  for 
bounding  the  channelized  receiver  detection  performance  of  the  methods  proposed 
in  Chapter  3.  The  desired  performance  for  any  robust  non-cooperative  technique 
would  be  better  than  that  of  a  radiometer,  with  somewhat  poorer  results  expected 
when  compared  with  matched  filter  detection  (assuming  AWGN  channel  conditions). 
Matched  filter  and  radiometric  probability  of  detection  (Pd)  results  for  PFA  =  10~2 
are  shown  in  Fig.  4.1  where  the  received  UWB  signal  is  as  shown  previously  in 
Fig.  3.1. 

4.3  Multi- Aperture  Cross  Correlation  Receiver  Detection  Performance 

Multi-aperture  cross  correlation  receiver  results  are  presented  for  completeness 
as  representing  one  means  for  performing  non-cooperative  detection  [16].  Probability 
of  detection  results  were  calculated  for  PFA  —  10-2  using  the  received  UWB  signal 
shown  in  Fig.  3.1.  As  shown  in  Fig.  4.2,  the  multi-aperture  cross  correlation  receiver 
of  Fig.  2.5  provides  an  improvement  of  approximately  1.5  dB  over  basic  radiometric 
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Figure  4.1  Matched  Filter  and  Radiometric  Detection  Results  for  Pfa  —  10  2 

detection  at  Pd  =  0.9.  These  results  are  consistent  with  the  theoretical  performance 
improvement  noted  in  [16]. 

4.4  Channelized  Receiver  Detection  Performance 

As  stated  in  Chapter  3,  channelized  receiver  detection  is  implemented  by 
performing  threshold  detection  on  the  processed  channelized  data  matrices  (TTM, 
CTM,  SSM,  and  CSM)  detailed  in  Section  3.2.  The  following  sections  provide  detec¬ 
tion  results  for  each  of  these  data  matrices  for  Pfa  —  10-2  using  the  received  UWB 
signal  in  Fig.  3.1. 

4.4.1  Temporal-Temporal  Matrix  (TTM)  Detection.  Probability  of 
detection  (Pd)  results  are  shown  in  Fig.  4.3  for  a  channelized  receiver  using  M  =  20, 
250  MHz  channels  and  a  varying  number  of  IFFT  points  ( Nifft )  to  form  the  TTM.  As 
indicated  in  Fig.  4.3,  detection  performance  improves  as  Ntjft  increases  until  Nt j ft 
becomes  greater  than  the  number  of  receiver  channels.  This  conclusion  is  supported 
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Figure  4.2  Cross  Correlation  Receiver  (CorrRx)  Detection  Results  Relative  to 
Matched  Filter  and  Radiometer  for  Pfa  —  10-2 

by  TTM  detection  results  for  channel  bandwidths  of  1.0  GHz,  500  MHz,  and  100  MHz 
as  included  in  Appendix  A.  For  a  channelized  receiver  employing  250  MHz  channels, 
the  best  performance  is  achieved  when  A, f ft  >  32,  which  corresponds  to  Nifft  values 
for  which  the  IFFT  no  longer  truncates  the  input  data  (note  that  A); j ft  —  32  through 
Nifft  —  256  data  is  coincident  in  Fig.  4.3). 

The  effect  of  varying  channel  bandwidth  is  considered  next.  The  number  of 
IFFT  points  was  selected  so  that  channelized  receivers  employing  different  band- 
widths  all  achieve  the  best  performance.  In  this  case,  N%fft  —  64  was  chosen  since 
the  smallest  channel  bandwidth  considered  is  100  MHz  which  requires  M  —  50  chan¬ 
nels  to  span  the  2.5  to  7.5  GHz  frequency  range.  Probability  of  detection  results 
for  different  receiver  channel  bandwidths  and  Nifft  —  64  are  shown  in  Fig.  4.4.  As 
indicated,  the  plots  for  all  four  channelized  receiver  bandwidths  are  perfectly  coin¬ 
cident,  indicating  that  detection  performance  is  independent  of  channel  bandwidth 
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Figure  4.3  TTM  Detection  Performance  for  a  Channelized  Receiver  with  M  —  20, 
250  MHz  Channels  and  Varying  Number  of  IFFT  Points 

when  Nifft  is  properly  chosen.  The  channelized  receiver  detection  performance  is 
approximately  2.5  dB  better  than  the  radiometer  at  Pd  —  0.9. 

4.4.2  Cross  Temporal  Matrix  (CTM)  Detection.  Probability  of  de¬ 
tection  (Pd)  results  are  shown  in  Fig.  4.5  for  a  channelized  receiver  using  M  —  20, 
250  MHz  channels  and  a  varying  number  of  IFFT  points  ( Nifft )  used  to  form  the 
CTM  from  the  corresponding  TTM.  As  with  TTM  detection,  CTM  detection  per¬ 
formance  improves  as  Nifft  increases  until  Ntf jt  becomes  greater  than  the  number 
of  channels.  This  conclusion  is  supported  by  the  CTM  results  in  Fig.  4.5,  as  well  as, 
additional  CTM  detection  results  for  channel  bandwidths  of  1.0  GHz,  500  MHz,  and 
100  MHz  as  included  in  Appendix  A.  Best  performance  is  achieved  when  Nifft  >  32, 
which  corresponds  to  Ntf ft  values  for  which  the  IFFT  no  longer  truncates  the  input 
data  (note  that  Nt jft  —  32  through  Nljft  —  256  data  is  coincident  in  Fig.  4.5). 

Paralleling  the  TTM  characterization  proceess,  the  effect  of  varying  channel 
bandwidth  is  considered  next.  The  number  of  IFFT  points  was  selected  so  that  chan- 
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Figure  4.4  TTM  Detection  Performance  for  a  Channelized  Receiver  Using 
Nifft  =  64  and  Varying  Channel  Bandwidth 
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Figure  4.5  CTM  Detection  Performance  for  a  Channelized  Receiver  with  M  =  20, 
250  MHz  Channels  and  Varying  Number  of  IFFT  Points 
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Figure  4.6  CTM  Detection  Performance  for  a  Channelized  Receiver  Using 
Nifft  —  64  and  Varying  Channel  Bandwidth 

nelized  receivers  employing  different  bandwidths  all  achieve  the  best  performance. 
Again,  Nifft  —  64  was  used  since  the  smallest  channel  bandwidth  considered  is 
100  MHz  which  requires  M  —  50  channels  to  span  the  2.5  to  7.5  GHz  frequency  range. 
Probability  of  detection  results  for  different  channelized  receiver  channel  bandwidths 
and  —  64  are  shown  in  Fig.  4.6.  These  results  indicate  that  increasing  receiver 
channel  bandwidth  improves  detection  performance.  The  performance  improvement 
over  the  radiometer  varies  from  approximately  2.0  dB  when  using  1.0  GHz  channels 
to  approximately  0.25  dB  when  using  100  MHz  channels  at  Pd  =  0.9. . 

4.4.3  Spectral-Spectral  Matrix  (SSM)  Detection.  Probability  of 
detection  (Pd)  results  are  shown  in  Fig.  4.7  for  a  channelized  receiver  with  M  —  20, 
250  MHz  channels  and  a  varying  number  of  FFT  points  ( Nfft )  used  to  form  the 
SSM.  Detection  performance  remains  approximately  the  same  for  all  Nfft  values 
considered.  This  conclusion  is  supported  by  SSM  results  in  Fig.  4.7,  as  well  as, 
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Figure  4.7  SSM  Detection  Performance  for  a  Channelized  Receiver  with  M  =  20, 
250  MHz  Channels  and  Varying  Number  of  FFT  Points 

additional  SSM  results  for  channel  bandwidths  of  1.0  GHz,  500  MHz,  and  100  MHz 
as  included  in  Appendix  A. 

The  effect  of  varying  the  channel  bandwidth  is  considered  next.  The  number 
of  FFT  points  was  selected  to  provide  both  reliable  detection  performance  and  ac¬ 
ceptable  simulation  run  times.  In  this  case,  Nfft  —  512  provided  adequate  results  in 
terms  of  both  performance  and  simulation  run  time  and  also  ensured  no  time  sam¬ 
ples  were  discarded  when  performing  the  FFT.  Probability  of  detection  results  for 
different  channel  bandwidths  and  Nfft  —  512  are  shown  in  Fig.  4.8.  Detection  perfor¬ 
mance  varies  slightly,  but  remains  approximately  the  same  over  the  range  of  channel 
bandwidths  considered.  Note  also  that  SSM  detection  performance  is  consistently 
poorer  than  radiometric  detection. 

4.4.4  Cross  Spectral  Matrix  (CSM)  Detection.  Probability  of  de¬ 
tection  (Pd)  results  are  shown  in  Fig.  4.9  for  a  channelized  receiver  using  M  —  20, 
250  MHz  channels  and  a  varying  number  of  FFT  points  ( Nfft )  used  to  form  the 
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Figure  4.8  SSM  Detection  Performance  for  a  Channelized  Receiver  Using 
N f ft  —  512  and  Varying  Channel  Bandwidth 

CSM  from  the  corresponding  SSM.  Detection  performance  remains  approximately 
the  same  for  all  N f ft  values  considered.  This  conclusion  is  supported  by  the  CSM 
results  in  Fig.  4.9,  as  well  as,  additional  CSM  results  for  channel  bandwidths  of 
1.0  GHz,  500  MHz,  and  100  MHz  as  included  in  Appendix  A. 

The  effect  of  varying  the  channel  bandwidth  is  considered  next.  The  number 
of  FFT  points  was  selected  to  provide  both  reliable  detection  performance  and  ac¬ 
ceptable  simulation  run  times.  In  this  case,  Nfft  —  512  provided  adequate  results 
in  terms  of  both  performance  and  simulation  run  time  and  also  ensured  no  time 
samples  were  discarded  when  performing  the  FFT.  Probability  of  detection  results 
for  different  channel  bandwidths  and  Nfft  —  512  are  shown  in  Fig.  4.10.  Unlike 
CTM  detection  performance  in  Fig.  4.6  which  shows  improvement  relative  to  the 
radiometer,  CSM  detection  performance  is  consistently  poorer  than  the  radiomter; 
CSM  performance  varies  slightly  but  remains  approximately  the  same  over  the  range 
of  bandwidths  considered. 
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SNR  (dB) 


Figure  4.9  CSM  Detection  Performance  for  a  Channelized  Receiver  with  M  —  20, 
250  MHz  Channels  and  Varying  Number  of  FFT  Points 


SNR  (dB) 

Figure  4.10  CSM  Detection  Performance  for  a  Channelized  Receiver  Using 
N f ft  —  512  and  Varying  Channel  Bandwidth 
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4.5  Downconverting  Channelized  Receiver  Detection  Performance 

As  stated  in  Chapter  3,  channelized  receiver  detection  is  implemented  by  per¬ 
forming  threshold  detection  on  the  processed  downconverting  channelized  receiver 
output  matrices.  Due  to  the  variation  in  detection  performance  of  the  downcon¬ 
verting  channelized  receiver  relative  to  the  downconversion  mixer  phases,  SNR  vs. 
Pd  results  for  the  channelized  receiver  processing  techniques  are  not  compared  with 
SNR  vs.  Pd  results  for  the  radiometer  and  the  matched  filter.  Instead,  the  detection 
performance  variation  is  plotted  as  a  function  of  the  initial  phases  on  the  downcon¬ 
version  mixers  for  a  constant  SNR  and  the  performance  of  the  proessing  technique 
is  compared  to  radiometric  performance  at  that  SNR.  All  detection  performance 
results  were  calculated  using  Pfa  —  10~2. 

4.5.1  Downconverting  Temporal-Temporal  Matrix  (TTM)  Detec¬ 
tion.  Probability  of  detection  (Pd)  results  are  shown  in  Fig.  4.11  for  a  downcon¬ 
verting  channelized  receiver  with  M  —  20,  250  MHz  channels  and  a  varying  number 
of  IFFT  points  ( Nifft )  used  to  form  the  TTM.  Detection  performance  improves  as 
Nifft  increases  and  appears  to  asymptotically  approach  a  limit  (i.e.  at  some  point 
increasing  Njjft  no  longer  improves  performance).  Considering  the  downcoverting 
TTM  results  in  Fig.  4.11,  along  with  additional  TTM  results  for  channel  bandwidths 
of  1.0  GHz,  500  MHz,  and  100  MHz  as  included  in  Appendix  A,  it  appears  that  best 
performance  is  achieved  when  Nifft  becomes  greater  than  twice  the  number  of  chan¬ 
nelized  receiver  channels,  or  Nifft  >  2 M .  It  is  noted  that  performance  for  N,j jt  —  32 
is  close  to  the  performance  for  Nifft  >  64  and  that  Ntfft  >  32  corresponds  to  Art j ft 
values  for  which  the  IFFT  no  longer  truncates  the  input  data. 

Detection  performance  response  due  to  variation  in  initial  downconversion 
mixer  phase  value  is  considered  next.  The  initial  phase  values  of  all  mixers  were 
set  equal  at  the  start  of  the  observation  interval  and  subsequently  varied  to  generate 
results  presented  in  Figs.  4.12  and  4.13.  The  signal-to-noise  ratio  (SNR)  was  held 
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Figure  4.11  TTM  Detection  Performance  for  a  Downconverting  Channelized  Re¬ 
ceiver  with  M  —  20,  250  MHz  Channels  and  Varying  Number  of  IFFT 
Points 

constant  at  SNR  —  0  dB  for  these  results.  As  illustrated  in  Fig.  4.13,  variation  in 
detection  performance  resulting  from  initial  phase  variation  decreases  as  the  number 
of  channels  is  increased  (corresponding  to  a  decrease  in  channelized  receiver  channel 
bandwidth). 

Radiometric  performance  at  SNR  —  0  dB  is  approximately  Pd  —  0.9  (see 
Fig.  4.1).  Thus,  data  in  Fig.  4.13  suggests  that  detection  using  TTM  data  can 
outperform  the  radiometer  when  using  1.0  GHz  and  500  MHz  channel  bandwidths. 
Detection  performance  varies  minimally  around  radiometric  performance  (Pd  =  0.9) 
for  channel  bandwidths  of  250  MHz  and  100  MHz. 

4.5.2  Downconverting  Cross  Temporal  Matrix  (CTM)  Detection. 

Probability  of  detection  (Pd)  results  are  shown  in  Fig.  4.14  for  a  downconverting 
channelized  receiver  with  M  —  20,  250  MHz  channels  and  a  varying  number  of 
IFFT  points  ( Nifft )  used  to  form  the  CTM  from  the  corresponding  TTM.  Detection 
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Figure  4.12 


Figure  4.13 


TTM  Detection  Performance  Versus  Mixer  Phase  for  a  Downconvert¬ 
ing  Channelized  Receiver  Using  SNR  —  0  dB  with  M  =  20,  250  MHz 
Channels 


TTM  Detection  Performance  Versus  Mixer  Phase  for  a  Downconvert¬ 
ing  Channelized  Receiver  Using  SNR  =  0  dB  with  Varying  Channel 
Bandwidth 
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Figure  4.14  CTM  Detection  Performance  for  a  Downconverting  Channelized  Re¬ 
ceiver  with  M  —  20,  250  MHz  Channels  and  Varying  Number  of  IFFT 
Points 

performance  improves  as  Nt f  ft  increases  until  Nt j ft  becomes  greater  than  the  number 
of  channels  used  in  the  channelized  receiver.  This  conclusion  is  supported  by  the 
downconverting  CTM  results  in  Fig.  4.14,  as  well  as,  additional  downconverting  CTM 
results  for  channel  bandwidths  of  1.0  GHz,  500  MHz,  and  100  MHz  as  included  in 
Appendix  A.  Best  performance  is  achieved  when  Nifft  >  32,  which  corresponds 
to  Nifft  values  for  which  the  IFFT  no  longer  truncates  the  input  data  (note  that 
Nifft  —  32  through  Nifft  —  256  data  is  coincident  in  the  figure). 

Detection  performance  response  due  to  variation  in  initial  downconversion 
mixer  phase  value  is  considered  next.  The  initial  phase  values  of  all  mixers  were 
set  equal  at  the  start  of  the  observation  interval  and  subsequently  varied  to  gener¬ 
ate  results  presented  in  Figs.  4.15  and  4.16.  The  signal-to-noise  ratio  (SNR)  was 
held  constant  at  SNR  =  0  dB  for  these  results.  Consistent  with  TTM  results  in 
Fig.  4.13,  CTM  results  in  Fig.  4.16  show  that  variation  in  detection  performance  re- 
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Figure  4.15  CTM  Detection  Performance  Versus  Mixer  Phase  for  a  Downconvert¬ 
ing  Channelized  Receiver  Using  SNR  =  0  dB  with  M  —  20,  250  MHz 
Channels 

suiting  from  initial  phase  variation  decreases  as  the  number  of  channels  is  increased 
(corresponding  to  a  decrease  in  channelized  receiver  channel  bandwidth). 

Radiometric  performance  at  SNR  —  0  dB  is  approximately  Pd  —  0.9  (see 
Fig.  4.1).  Thus,  data  in  Fig.  4.16  suggests  that  detection  using  CTM  data  can 
outperform  the  radiometer  when  using  a  channel  bandwidth  of  1.0  GHz.  Detection 
performance  exceeds  radiometric  performance  and  shows  very  little  variation  due  to 
phase  for  channel  bandwidths  of  500  MHz,  250  MHz,  and  100  MHz. 

4.5.3  Downconverting  Spectral-Spectral  Matrix  (SSM)  Detection. 

Probability  of  detection  (Pd)  results  are  shown  in  Fig.  4.17  for  a  downconvert¬ 
ing  channelized  receiver  with  M  —  20,  250  MHz  channels  and  a  varying  number  of 
FFT  points  ( Nfft )  used  to  form  the  SSM.  Detection  performance  remains  approxi¬ 
mately  the  same  for  all  Nfft  values  considered.  This  conclusion  is  supported  by  the 
downconverting  SSM  results  in  Fig.  4.17,  as  well  as,  additional  downconverting  SSM 
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Figure  4.16  CTM  Detection  Performance  Versus  Mixer  Phase  for  a  Downconvert¬ 
ing  Channelized  Receiver  Using  SNR  =  0  dB  with  Varying  Channel 
Bandwidth 

results  for  channel  bandwidths  of  1.0  GHz,  500  MHz,  and  100  MHz  as  included  in 
Appendix  A. 

Detection  performance  response  due  to  variation  in  initial  downconversion 
mixer  phase  value  is  considered  next.  The  initial  phase  values  of  all  mixers  were 
set  equal  at  the  start  of  the  observation  interval  and  subsequently  varied  to  generate 
results  presented  in  Figs.  4.18  and  4.19.  The  SNR  was  kept  constant,  at  SNR  —  0  dB 
and  the  number  of  FFT  points  was  chosen  to  be  Nf  jt  =  512  (to  ensure  no  time  sam¬ 
ples  are  discarded  when  performing  the  FFT).  As  illustrated  in  Fig.  4.19,  variation 
in  detection  performance  resulting  from  phase  variations  decreases  slightly  as  the 
number  of  channels  increases.  However,  the  detection  performance  variation  is  not 
reduced  to  that  exhibited  in  the  downconverting  TTM  and  CTM  results  of  Fig.  4.11 
and  Fig.  4.14,  respectively. 
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Figure  4.17 


Figure  4.18 


SSM  Detection  Performance  for  a  Downconverting  Channelized  Re¬ 
ceiver  with  M  —  20,  250  MHz  Channels  and  Varying  Number  of  FFT 
Points 


SSM  Detection  Performance  Versus  Mixer  Phase  for  a  Downconverting 
Channelized  Receiver  Using  SNR  =  0  dB  with  M  —  20,  250  MHz 
Channels 
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Figure  4.19  SSM  Detection  Performance  Versus  Mixer  Phase  for  a  Downconvert¬ 
ing  Channelized  Receiver  Using  SNR  —  0  dB  with  Varying  Channel 
Bandwidth 

Given  radiometric  performance  at  SNR  —  0  dB  is  approximately  Pd  =  0.9, 
results  in  Fig.  4.19  clearly  indicate  that  downconverting  SSM  detection  performance 
is  much  poorer  than  the  radiometer  for  all  channel  band  widths  considered. 

4.5.4  Downconverting  Cross  Spectral  Matrix  (CSM)  Detection. 

Probability  of  detection  (Pd)  results  are  shown  in  Fig.  4.20  for  a  downconvert¬ 
ing  channelized  receiver  with  M  —  20,  250  MHz  channels  and  a  varying  number 
of  FFT  points  ( Nfjt )  used  to  form  the  CSM  from  the  corresponding  SSM.  Detec¬ 
tion  performance  remains  approximately  the  same  for  all  Nfft  values  considered. 
This  conclusion  is  supported  by  the  downconverting  CSM  results  in  Fig.  4.20,  as 
well  as,  additional  downconverting  CSM  results  for  channel  bandwidths  of  1.0  GHz, 
500  MHz,  and  100  MHz  as  included  in  Appendix  A. 

Detection  performance  response  due  to  variation  in  initial  downconversion 
mixer  phase  value  is  considered  next.  The  initial  phase  values  of  all  mixers  were  set 
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Figure  4.20  CSM  Detection  Performance  for  a  Downconverting  Channelized  Re¬ 
ceiver  with  M  —  20,  250  MHz  Channels  and  Varying  Number  of  FFT 
Points 

equal  at  the  start  of  the  observation  interval  and  subsequently  varied  to  generate  re¬ 
sults  presented  in  Figs.  4.21  and  4.22.  The  SNR  was  kept  constant,  at  SNR  =  0  dB 
and  the  number  of  FFT  points  was  chosen  to  be  Nfft  —  512  (to  ensure  no  time 
samples  are  discarded  when  performing  the  FFT).  As  illustrated  in  Fig.  4.22,  the 
variation  in  detection  performance  due  to  phase  variations  decreases  as  the  number 
of  channels  is  increased,  and  at  the  same  time  detection  performance  increases. 

Given  radiometric  performance  at  SNR  —  0  dB  is  approximately  Pd  =  0.9, 
data  in  Fig.  4.22  suggests  that  downconverting  CSM  detection  outperforms  the 
radiometer  when  using  a  receiver  channel  bandwidth  of  100  MHz.  Radiometric 
detection  outperforms  CSM  detection  for  all  other  channel  bandwidths  considered 
(1.0  GHz,  500  MHz,  and  250  MHz)  with  performance  improving  as  channel  band¬ 
width  decreases  (corresponding  to  an  increase  in  the  number  of  channels). 
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Figure  4.21 


Figure  4.22 


Phase  (deg) 


CSM  Detection  Performance  Versus  Mixer  Phase  for  a  Downconvert¬ 
ing  Channelized  Receiver  Using  SNR  —  0  dB  with  M  —  20,  250  MHz 
Channels 


CSM  Detection  Performance  Versus  Mixer  Phase  for  a  Downconvert¬ 
ing  Channelized  Receiver  Using  SNR  =  0  dB  with  Varying  Channel 
Bandwidth 
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4.6  Summary 

Detection  performance  results  were  provided  for  each  of  the  receivers  and 
processing  techniques  introduced  in  Chapters  2  and  3.  First,  matched  filter  and 
radiometric  detection  performance  results  were  introduced  for  the  UWB  waveform 
under  consideration.  These  results  effectively  bound  anticipated  results  for  the  pro¬ 
posed  non-cooperative  techniques.  For  completeness,  detection  performance  of  a 
multi-aperture  cross  correlation  technique  was  introduced  next  as  being  representa¬ 
tive  of  one  specific  non-cooperative  detection  technique.  Finally,  detection  perfor¬ 
mance  results  using  the  matrix  data  (TTM,  CTM,  SSM  and  CSM)  and  techniques 
introduced  in  Chapter  3  were  presented  for  both  the  channelized  receiver  and  the 
downconverting  channelized  receiver. 
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5.  Conclusions 


5.1  Summary 

The  majority  of  this  thesis  has  been  devoted  to  introducing  and  analyzing 
channelized  receiver  processing  techniques  for  use  in  detecting  ultra  wideband  sig¬ 
nals,  specifically,  a  single  UWB  impulse  in  noise.  The  detection  techniques  considered 
operated  on  data  contained  in  the  temporal-temporal  matrix  (TTM),  created  by  per¬ 
forming  an  IFFT  on  columns  of  the  M  x  S  Channelized  Data  Matrix  (CDM);  the 
cross  temporal  matrix  (CTM),  created  by  performing  a  correlation  between  all  com¬ 
binations  of  TTM  columns;  the  spectral-spectral  matrix  (SSM),  created  by  perform¬ 
ing  an  FFT  on  the  rows  of  the  CDM;  and  the  cross  spectral  matrix  (CSM),  created 
by  performing  a  correlation  between  all  combinations  of  SSM  columns.  Threshold 
detection  was  performed  using  each  of  these  matrices  and  probability  of  detection 
results  generated  for  varying  parameters,  including  received  SNR,  the  number  of 
IFFT  and  FFT  points  used  to  form  the  TTM  (CTM)  and  SSM  (CSM),  respectively, 
receiver  channel  bandwidth,  and  the  initial  phase  value  of  downconversion  mixers 
(downconverting  channelized  receiver  only). 

Three  other  receiver  architectures  were  considered,  including  the  matched  fil¬ 
ter,  the  radiometer,  and  the  multi-aperture  cross  correlation  receiver.  Detection 
performance  of  the  matched  filter  and  radiometric  receivers  serve  as  performance 
bounds  (upper  and  lower,  respectively)  for  characterizing  channelized  receiver  per¬ 
formance.  Detection  performance  of  the  multi-aperture  cross  correlation  receiver  is 
provided  by  way  of  completeness  in  that  its  performance  is  representative  of  what 
can  be  achieved  using  non-cooperative  (non-matched  filter)  techniques. 

5.2  Conclusions 

5.2.1  Channelized  Receiver  Detection  Performance.  Detection  per¬ 
formance  of  the  channelized  receiver  using  the  TTM  (with  an  optimum  number  of 
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IFFT  points  as  discussed  in  Chapter  4)  provides  approximately  2.5  dB  improvement 
over  the  radiometer  while  performing  approximately  3.5  dB  poorer  than  the  matched 
filter  at  Pd  —  0.9.  Likewise,  detection  performance  using  the  CTM  is  better  than 
radiometric  performance  (using  the  optimum  number  IFFT  points)  but  not  quite  as 
good  as  that  realized  when  using  the  TTM.  It  was  also  shown  that  CTM  detection 
performance  degrades  as  channel  bandwidth  decreases. 

Detection  performance  using  either  of  the  spectral-based  processing  techniques 
(i.e.  SSM  or  CSM)  was  poorer  than  that  of  the  radiometer.  The  best  performance 
achieved  by  either  of  these  techniques  was  approximately  4.0  dB  poorer  than  radio- 
metric  detection  at  Pd  —  0.9.  This  result  was  expected  since  the  power  in  a  UWB 
signal  is  spread  across  an  extremely  wide  bandwidth,  as  shown  in  Fig.  2.2,  such  that 
channel  noise  quickly  overcomes  the  UWB  signal  in  the  spectral  domain  as  the  SNR 
decreases. 

5.2.2  Downconverting  Channelized  Receiver  Performance.  Detec¬ 
tion  performance  of  the  downconverting  channelized  receiver  using  any  of  the  chan¬ 
nelized  receiver  processing  techniques  varies  relative  to  the  downconversion  mixer 
phases.  For  detection  using  the  TTM,  detection  performance  variation  decreases 
as  the  number  of  channels  increases  (channel  bandwidth  decreases)  and  for  chan¬ 
nel  bandwidths  of  250  MHz  and  100  MHz,  detection  performance  approximately 
equals  that  of  the  radiometer.  Detection  performance  using  the  CTM  generated 
with  1.0  GHz  channels  provides  detection  ranging  from  approximately  Pd  —  0.7  to 
Pd  =  1  as  mixer  phase  varied  from  0  to  360  degrees.  For  downconverting  channelized 
receivers  using  100  MHz,  250  MHz,  or  500  MHz  channels,  detection  performance  sur¬ 
passes  that  of  the  radiometer  while  remaining  relatively  constant  for  all  phase  values. 

Variation  in  downconversion  mixer  phase  resulted  in  SSM  detection  perfor¬ 
mance  variation  for  all  channel  bandwidths,  although  the  amount  of  variation  in  de¬ 
tection  performance  appeared  to  decrease  as  the  number  of  channels  was  increased 
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(decreasing  channel  bandwidth).  Detection  performance  using  the  SSM  does  not 
come  close  to  radiometric  performance.  Using  the  CSM,  detection  performance  vari¬ 
ation  as  a  function  of  mixer  phase  occurs  which  decreases  as  the  number  of  channels 
increases.  The  CSM  detection  performance  surpasses  that  of  the  radiometer  for 
the  downconverting  channelized  receiver  using  100  MHz  channels,  in  which  case  the 
performance  variation  as  a  function  of  mixer  phase  was  minimal. 

5.3  Recommendations  for  Future  Research 

5.3.1  Detection  Performance  in  Coexistence  Scenarios.  This  work 
only  considered  received  signals  consisting  of  a  single  UWB  signal  with  and  without 
additive  white  Gaussian  channel  noise.  Unfortunately,  real  world  environments  can 
contain  many  different  signals  using  many  different  modulation  types.  Performance 
of  the  channelized  receiver  processing  techniques  presented  here  could  be  analyzed 
in  coexisting  scenarios,  i.e. ,  consider  received  signals  containing  the  UWB  signal 
considered  here  along  with  other  signals  such  as  narrowband  tones,  M-ary  phase 
shift  keyed  (MPSK),  M-ary  frequency  shift  keyed  (MFSK),  direct  sequence  spread 
spectrum  (DSSS),  and  other  modulation  types  encountered  in  communication,  nav¬ 
igation,  and  radar  applications. 

5.3.2  Channel  Assessment  and  Characterization.  In  some  cases, 
it  may  be  beneficial  to  identify  the  number,  type,  and  nature  of  signals  present  in 
an  environment  at  any  given  time  (channel  assessment).  With  narrowband  (char¬ 
acterized  relative  to  UWB)  signals  present  in  the  environment,  it  may  be  difficult 
to  reliably  detect  and/or  characterize  UWB  signals.  Visual  inspection  of  the  two 
dimensional  TTM,  CTM,  SSM,  and  CSM  plots  in  Chapter  3  suggest  the  possibil¬ 
ity  of  employing  pattern  recognition  techniques  in  conjunction  with  the  channelized 
receiver  processing  techniques  presented  in  this  work  to  make  this  determination. 
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Appendix  A. 


A.l  Channelized  Receiver  Detection  Performance 


Figure  A.l  TTM  -  1  GHz  ChRx  (Varying  Number  of  Points  in  IFFT) 
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A. 2  Down-Converting  Channelized  Receiver  Detection  Performance 
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Figure  A.  13  TTM  -  1  GHz  DC  ChRx  (Varying  Number  of  Points  in  IFFT) 
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Appendix  B.  Alternate  UWB  Pulse  Offset 


B.l  Input  Signal  and  Processed  Matrices 


Figure  B.l  UWB  ChRx  Input  Signal 
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Figure  B.ll  TTM  -  Nifft  =  64 
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Figure  B.12  CTM  -  250  MHz  ChRx  (Varying  Number  of  Points  in  IFFT) 
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Figure  B.13  CTM  -  Nt j ft  =  64 
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Figure  B.14  SSM  -  250  MHz  ChRx  (Varying  Number  of  Points  in  IFFT) 


B.3  Downconverting  Channelized  Receiver  Detection  Performance 
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Figure  B.18 
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TTM  -  250  MHz  DC  CliRx  (Varying  Number  of  Points  in  IFFT) 
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Figure  B.19 
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TTM  -  250  MHz  DC  ChRx,  SNR  =  0  dB  (Varying  Starting  Phases) 


Figure  B.20  TTM  -  DC  ChRx,  SNR  =  0  dB  (Varying  Starting  Phases) 
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Figure  B.23  CTM  -  DC  ChRx,  SNR  =  0  dB  (Varying  Starting  Phases) 
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Figure  B.24  SSM  -  250  MHz  DC  ChRx  (Varying  Number  of  Points  in  IFFT) 
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Figure  B.25 
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SSM  -  250  MHz  DC  ChRx,  SNR  =  0  dB  (Varying  Starting  Phases) 


Figure  B.26  SSM  -  DC  ChRx,  SNR  =  0  dB  (Varying  Starting  Phases) 


Appendix  C.  MATLAB  Code 
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7o  UWB  Power  Spectral  Density 
7.  —  Analytic 

0  /  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

clear  all,  close  all,  clc,  format  compact 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o 


7o  UWB  Signal  Parameters 

07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


fc  =  5e9; 

Tw  =  2/fc; 

Tm  =  0.4*Tw; 
Ts  =  2*Tw; 
delta  =  Tw/2; 


7.  center  frequency 
7o  pulse  duration 
7.  pulse  width  parameter 
7o  symbol  duration 
7o  position  modulation  offset 


f  =  0 : 1/ (1000*Tw) : 40e9 ;  %  frequency  vector 


7o  FT  of  received  UWB  waveform 

Wf  =  pi*Tm~3*f  .  '"2/sqrt  (2)  .  *exp(-pi/2*Tm~2*f  .  ~2)  ; 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 
/o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7o  Analytic  PSD  —  Uniform 

07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0 J  0  7  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

Suni_con  =  le-13*ones (1 , length(f ) )  ;  °/„ 
Suni_dis  =  l/Ts~2*ones  (1 ,  length(f ) ) ;  °/„ 


n  =  f ind(mod(f , 1/Ts)  ~=  0); 
for  i  =  l:length(n), 

Suni_dis (n(i) )  =  le-13;  % 

end 

Smod_uni  =  Suni_con  +  Suni_dis;  % 

Suni  =  Smod_uni  .  *  (Wf  .  *conj  (Wf ) )  ;  °/„ 

Suni(l)  =  le-13;  % 

Suni_dB  =  10*logl0 (abs (Suni) ) ;  % 


continuous  part 
discrete  part 

’ discretize'  discrete  part 

sum  continuous  and  discrete  parts 
total  signal  PSD 
avoid  dB  conversion  warning 
total  signal  PSD  (dB) 


figure (1) 

plot (f , Suni_dB , J  kx ’ ) 
%axis ( [0  16e9  -160  -10]) 
ylabel(’ Magnitude  (dB) ’ ) 
xlabel( ’Frequency  (Hz)’) 


C-l 


0,  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  70  /o  7o  7o  7o  7o  7o  7o  7o  7o  7o  7o  7o  7o  7o  7o  7o  7o  7o 


%  Analytic  PSD  —  PPM 

07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 
70  70  /0  /0  /O  /O  70  /O  70  70  /O  /O  /O  /O  70  /O  70  70  /O  /O  /O  70  70  70  70 

Sppm_con  = 

Sppm_dis  = 
n  =  f ind(mod(f , 1/Ts)  ~=  0); 
for  i  =  l:length(n), 

Sppm_dis (n(i) )  =  le-13; 

end 

Smod_ppm  =  Sppm_con  +  Sppm_dis; 
Sppm  =  Smod_ppm  . *  (Wf . *conj (Wf ) ) ; 
Sppm(l)  =  le-13; 

Sppm_dB  =  10*logl0 (abs (Sppm) ) ; 


continuous  part 
discrete  part 

%  'discretize'  discrete  part 

%  sum  continuous  and  discrete  parts 
7o  total  signal  PSD 
7o  avoid  dB  conversion  warning 
7.  total  signal  PSD  (dB) 


(1/ (2*Ts) ) * (1-cos (4*pi*f *delta) ) ;  % 

(1/ (2*Ts~2) ) * (1+cos (4*pi*f *delta) ) ;  % 


figure (1) ,hold  on 
plot (f , Sppm_dB , ' k ' ) 
axis ( [0  18e9  -180  -10]); 
ylabel ( 'Magnitude  (dB)’); 
xlabel (’ Frequency  (Hz)’); 
grid 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7.  Analytic  PSD  —  PAM 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


Spam_con  =  2/Ts; 

Spam_dis  =  0; 

Smod_pam  =  Spam_con  +  Spam_dis; 
Spam  =  Smod_pam  . *  (Wf . *conj (Wf ) ) ; 
Spam(l)  =  le-13; 

Spam_dB  =  10*logl0 (abs (Spam) ) ; 


7o  continuous  part 
7.  discrete  part 

7.  sum  continuous  and  discrete  parts 
7.  total  signal  PSD 
7.  avoid  dB  conversion  warning 
7.  total  signal  PSD  (dB) 


figure (1) 

plot (f , Spam_dB , ' k — ' ) 

7.  axis  ( [0  16e9  -160  -100]); 
ylabel (’Magnitude  (dB) ’ ) ; 
xlabel ( ’ Frequency  (Hz) ’ ) ; 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  Analytic  PSD  —  BPPM 

07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


Sbppm_con  =  1/Ts; 


7o  continuous  part 
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Sbppm_dis  =0;  %  discrete  part 

Smod_bppm  =  Sbppm_con  +  Sbppm_dis;  sum  continuous  and  discrete  parts 
Sbppm  =  Smod_bppm  .  *  (Wf . *conj (Wf)) ;%  total  signal  PSD 
Sbppm(l)  =  le-13;  %  avoid  dB  conversion  warning 

Sbppm_dB  =  10*logl0 (abs (Sbppm) ) ;  total  signal  PSD  (dB) 

figure (1) 

plot (f , Sbppm_dB , ’ k : ’ ) 

*/.  axis  ( [0  16e9  -160  -100]); 
ylabel(’ Magnitude  (dB) ’ ) ; 
xlabel ( ’ Frequency  (Hz) ’ ) ; 

legend ( ’ Uniform ‘ , 'PPM1 , 'PAM' , 'BPPM') 
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%  Brett  D .  Gronholz 
%  EENG  799  —  Sunnner/Fall  2004 
7o  UWB  Power  Spectral  Density 
7o  —  Simulated 

y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

clear  all,  close  all,  clc,  format  compact 


»/  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  UWB  Signal  Parameters 

07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  / 0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  / 0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


fc  =  5e9; 

Tw  =  2/fc; 

Ts  =  2*Tw; 

To  =  Ts; 
delt  =  0.04*Tw; 
Ns  =  1000; 

P  =  1; 
jtr  =  0; 
dly  =  0; 


7.  center  frequency 

7,  pulse  duration 

7,  symbol  duration 

7.  symbol  repetition  interval 

7o  time  resolution 

7.  number  of  symbols 

7.  signal  power 

7o  jitter  as  percentage  of  Ts 

7.  first  pulse  delay 


07070707070707070707070707 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7o  Uniform 

07070707070707070707070707 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


uni  =  uwb(Tw,To,delt,Ns,P, jtr, ’uni 
%  uni  =  awgn(uni , -10 , ’measured’ ) ; 


>  dly)  ; 


%  uniform  pulse  train 
7,  add  AWGN 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o 


7o  PSD  of  uniform  pulse  train 

07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


uni_tot  =  length(uni) ; 
f  =  l/delt*(0:uni_tot/2)/uni_tot; 
UNI  =  fft(uni); 

Suni  =  (abs (UNI/uni_tot) . ~2) *Ns*To ; 
for  i  =  1 : length(Suni) , 
if  Suni(i)  ==  0, 

Suni(i)  =  le-13; 

end 


end 

Suni_dB  =  10*logl0 (Suni) ; 


7,  length  of  signal 
7o  frequency  vector 
7o  FFT 
7o  PSD 


7o  avoid  dB  conversion  warning 

7,  PSD  (dB) 


figure (1) 

plot (f ,Suni_dB(l :uni_tot/2+l) , ’bx’ ) 
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axis ( [0  20e9  -170  -40]); 
ylabel ( ’Magnitude  (dB)’); 
xlabel (’ Frequency  (Hz)’); 
grid 


yyyyyyyyy 

/o  /o  /o  /o  /o  /o  /o  /o  /o 

7a  PPM 
yyyyyyyyy 

/o  /o  /o  /o  /o  /o  /  0  /o  /o 

ppm  =  uwb(Tw,To,delt,Ns,P,  jtr,  ’ppm’  ,dly)  ;  7„  PPM  pulse  train 

7a  ppm  =  awgn(ppm, -10 ,  ’measured’ )  ;  7»  add  AWGN 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7a  PSD  of  PPM  Pulse  Train 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o 


ppm_tot  =  length (ppm) ; 
f  =  l/delt*(0:ppm_tot/2)/ppm_tot; 
PPM  =  fft(ppm); 

Sppm  =  (abs (PPM/ppm_tot) . ~2) *Ns*To ; 
Sppm_dB  =  10*logl0 (Sppm) ; 


7a 

7. 

7a 

7. 

7a 


length  of  signal 
frequency  vector 
FFT 
PSD 

PSD  (dB) 


figure (2) 

plot (f ,Sppm_dB(l :ppm_tot/2+l) , ’b ’ ) 
axis ( [0  20e9  -170  -40]); 
ylabel ( ’ Magnitude  (dB) ’ ) ; 
xlabel (’ Frequency  (Hz)’); 
grid 


yyyyyyyyy 

/o  /o  /o  /o  /o  /o  /o  /o  /o 

7a  PAM 


yyyyyyyyy 

/o  /o  /o  /o  /o  /o  /o  /o  /o 


pam  =  uwb(Tw,To,delt,Ns,P, jtr, ’pam’ 
7a  pam  =  awgn(pam, -10 , ’measured’ )  ; 


>  dly)  ; 


7.  PAM  pulse  train 
7o  add  AWGN 


■ /  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7a  PSD  of  PAM  Pulse  Train 
o/  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

pam_tot  =  length (pam) ; 
f  =  l/delt*(0:pam_tot/2)/pam_tot; 

PAM  =  fft(pam); 

Spam  =  (abs (PAM/pam_tot) . ~2) *Ns*2*To ; 
Spam_dB  =  10*logl0 (Spam) ; 


7. 

7. 

7. 

7a 

7. 


length  of 
frequency 
FFT 
PSD 

PSD  (dB) 


signal 

vector 
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figure (3) 

plot (f ,Spam_dB(l :pam_tot/2+l) , ’b ’ ) 
axis ( [0  20e9  -170  -40]); 
ylabel ( ’Magnitude  (dB)’); 
xlabel( ’Frequency  (Hz)’); 
grid 


yyyyyyyyyy 

/o  /o  /o  /o  /  0  /o  /o  /o  /o  /o 

%  BPPM 
yyyyyyyyyy 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

bppm  =  uwb(Tw,To,delt,Ns,P, jtr, ’bppm’ ,dly) 
%  bppm  =  awgn (bppm, -10, ’measured’) ; 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


°/a  PSD  of  BPPM  pulse  train 
y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


bppm_tot  =  length (bppm) ; 
f  =  l/delt*(0:bppm_tot/2)/bppm_tot; 
BPPM  =  fft(bppm); 

Sbppm  =(abs(BPPM/bppm_tot) . ~2) *Ns*To ; 
Sbppm_dB  =  10*logl0 (Sbppm) ; 


7. 

7o 

7. 

1 

1 


figure (4) 

plot (f ,Sbppm_dB(l :bppm_tot/2+l) , ’b’ ) 

axis ( [0  20e9  -170  -40]); 

ylabel (’Magnitude  (dB) ’ ) ; 

xlabel(’ Frequency  (Hz)’); 

grid 


;  7o  BPPM  pulse  train 
7o  add  AWGN 


length  of  signal 
frequency  vector 
FFT 
PSD 

PSD  (dB) 
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%  Brett  D .  Gronholz 
%  EENG  799  -  Summer/Fall  2004 
%  Wideband  Correlation  Receiver 
%  -  UWB  Signal  Detection 

y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o 

clear  all  close  all,  clc,  format  compact 


»/  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  Simulation  Parameters 


07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o 


SNR  =  0; 
fl  =  2 . 5e9 ; 
fh  =  7 . 5e9 ; 
Nz  =  400; 

N  =  4; 


7.  input  signal  SNR  in  dB 

7,  lower  BPF  frequency 

7o  upper  BPF  frequency 

7.  one-sided  zero  padding  length  for 

7„  filter  order 


’f iltfilt’ 


07  07  07  07  0  7  0  7  07  07  07  07  07  07  0  7  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  UWB  Signal  Parameters 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


fc  =  5e9; 

Tw  =  2/fc; 

Ts  =  2*Tw; 

To  =  10*Ts; 
delt  =  0. Ole-9; 
fs  =  1/delt; 

Ns  =  1; 

P  =  1; 
jtr  =  0; 
method  =  'uni'; 
dlyl  =  2e-9; 
dly2  =  4e-9; 


7,  center  frequency 

7o  pulse  duration 

7.  symbol  duration 

7.  symbol  repetition  interval 

7o  time  resolution 

7.  sample  frequency 

7.  number  of  symbols 

7.  signal  power 

7.  jitter  as  percentage  of  Ts 

7o  modulation  method 

7.  first  pulse  delay,  sig  1 

7.  first  pulse  delay,  sig  2 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7.  Generate  Signals 

07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o 


xl  =  uwb(Tw, To, delt, Ns, P, jtr, method, dlyl) ; 
x2  =  uwb(Tw, To, delt, Ns, P, jtr, method, dly2) ; 
Pxl  =  sum(xl . ~2) /length(xl) ; 

Px2  =  sum(x2 . ~2) /length(x2) ; 
nl  =  randn(l , length(xl) ) ; 
n2  =  randn(l , length (x2) ) ; 


7o  UWB  signal  out  of  antenna  1 
7o  UWB  signal  out  of  antenna  2 
7o  power  in  UWB  sig  1 
7o  power  in  UWB  sig  2 
7.  noise  realization  1 
7.  noise  realization  2 
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y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o/o  /o  /o  /o 

%  Filter  Input  Signals 
■ /  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o 

[b,a]  =  butter(N, [f l/(f s/2)  fh/(f s/2)] ) ; % 


xfl  =  [zeros(l,Nz)  xl  zeros(l,Nz)] 
xfl  =  real(f iltf ilt (b, a,xf 1)) ; 
xfl  =  xf 1 (Nz+1 : end-Nz) ; 

Pxfl  =  sum(xf 1 . "2) /length(xf 1) ; 

xf2  =  [zeros(l,Nz)  x2  zeros(l,Nz)] ; 
xf2  =  real(f iltfilt(b,a,xf2)) ; 
xf2  =  xf 2 (Nz+1 : end-Nz) ; 

Pxf2  =  sum(xf 2 . "2) /length(xf 2) ; 

nfl  =  [zeros(l,Nz)  nl  zeros(l,Nz)] ; 

nfl  =  filtf ilt(b,a,nf 1) ; 

nfl  =  nfl (Nz+1 : end-Nz) ; 

nfl  =  nf 1/sqrt (var (nf 1)) ; 

nfl  =  sqrt (Pxfl/ 10~ (SNR/ 10) ) *nf 1 ; 

nf2  =  [zeros(l,Nz)  n2  zeros(l,Nz)] ; 

nf2  =  filtf ilt(b,a,nf2) ; 

nf2  =  nf 2 (Nz+1 : end-Nz) ; 

nf2  =  nf 2/sqrt (var (nf 2) ) ; 

nf 2  =  sqrt (Pxf 2/10~ (SNR/10) ) *nf 2 ; 


filter  coeffs 

zero-pad 
l  filter 
%  remove  zeros 

l  power  in  filtered  UWB  signal 

%  zero-pad 
%  filter 
l  remove  zeros 

/  power  in  filtered  UWB  signal 

%  zero-pad 
%  filter 
%  remove  zeros 
%  normalize  noise  power 
l  noise  at  required  SNR 

*/.  zero-pad 
%  filter 
%  remove  zeros 
70  normalize  noise  power 
70  noise  at  required  SNR 


7. 


0  /  0  /  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o 


7o  Correlation 


07  07  07  07  0  7  0  7  07  07  07  07  07  07  0  7  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


[xc,lags]  =  xcorr(xf 1+nf 1 ,xf2+nf2) ; 
plot (lags , abs (xc) ) 
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%  Brett  D .  Gronholz 
%  EENG  799  —  Sunnner/Fall  2004 
%  UWB  Detection  Probability  (Pd) 
%  -  Cross-Correlation  Receiver 


0/  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


clear  all,  close  all,  clc,  format  compact 


07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7o  Simulation  Parameters 

07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


Wc  =  [le9] ; 

7.  channel  bandwidths  to  simulate 

Nz  =  [500] ; 

7.  zero-padding  lengths  for  each  ’Wc’  (one-sided) 

fl  =  2 . 5e9 ; 

7.  lower  Hrf/ChRx  frequency 

fh  =  7 . 5e9 ; 

7.  upper  Hrf/ChRx  frequency 

o 

II 

u 

T3 

7.  downconvert?  (l=yes,0=no) 

ph  =  0; 

7.  DC  phases  to  simulate 

meth  =  1 ; 

7,  processing  method  (1=TTM,2=SSM,3=CTM,4=CSM) 

Nfft  =  64; 

7.  (I)FFT  length 

Pfa  =  10~-2 ; 

7.  probability  of  false  alarm 

R  =  10/Pfa; 

7.  number  of  realizations 

SNR  =  [-10:1:5] ; 

7.  SNRs  to  simulate 

o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  UWB  Signal  Parameters 


07  07  07  07  07  07  07  07  07  07  07  07  I J  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


fc  =  5e9; 

Tw  =  2/fc; 

Ts  =  2*Tw; 

To  =  Ts/2 ; 
delt  =  0. Ole-9 
fs  =  1/delt; 

Ns  =  1; 

P  =  1; 
jtr  =  0; 
method  =  'uni’ 
dlyl  =0;  l 

dly2  =  0;  7. 


7  center  frequency 
7.  pulse  duration 
7o  symbol  duration 
7o  symbol  repetition  interval 
7.  time  resolution 
7o  sample  frequency 
7o  number  of  symbols 
l  signal  power 

1  jitter  as  percentage  of  Ts 
7.  modulation  method 
first  pulse  delay,  sig  1 
first  pulse  delay,  sig  2 


07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7.  Generate  Signals 

07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0 


xl  =  uwb(Tw, To, delt, Ns, P, jtr, method, dlyl) ; 


7.  UWB  signal  at  output  of  antenna  1 
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xl  =  [zeros (1 , 240)  xl  zeros (1, 240) ] ; 

x2  =  uwb(Tw, To, delt, Ns, P,jtr, method, dly2) ; 

x2  =  [zeros(l, 100)  x2  zeros(l,380)]  ; 

Pxl  =  sum(xl . ~2) /length(xl) ; 

Px2  =  sum(x2 . ~2) /length(x2) ; 
nl  =  randn(R, length(xl) ) ; 
n2  =  randn(R, length (x2) ) ; 


7.  add  zeros 

7.  UWB  signal  at  output  of  antenna  2 

%  add  zeros 

7»  power  in  UWB  sig  1 

7»  power  in  UWB  sig  2 

7»  generate  noise  1 

7o  generate  noise  2 


»/  0  7  07  07  07  0 I  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  Filter  Input  Signals 


07  07  07  07  0  7  0  7  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  0  7  0  7  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


N  =  4;  7,  order  of  Hrf  BPF 

NzRF  =  200;  °/„  one-sided  zero  padding  length  for 

[b,a]  =  butter  (N,  [f  1/ (f  s/2)  fh/(fs/2)] ) filter  coeffs 


’filtfilt’ 


xfl  =  [zeros (1 , NzRF)  xl  zeros (1 , NzRF)] ; 
xfl  =  real(f iltf ilt (b, a,xf 1)) ; 
xfl  =  xf 1 (NzRF+1 : end-NzRF) ; 

Pxfl  =  sum(xf 1 . "2) /length(xf 1) ; 


7o  zero-pad 

%  filter 

%  remove  zeros 

*/.  power  in  filtered  UWB  sig 


xf2  =  [zeros (1 , NzRF)  x2  zeros (1 , NzRF)] ; 
xf2  =  real (filtfilt (b, a, xf 2) ) ; 
xf 2  =  xf 2 (NzRF+1 : end-NzRF) ; 

Pxf2  =  sum(xf 2 . ~2) /length(xf 2) ; 


7o  zero-pad 

7,  filter 

7.  remove  zeros 

7o  power  in  filtered  UWB  sig 


ntl  =  [zeros (R, NzRF)  nl  zeros (R, NzRF) ] ; 
for  i  =  1:R, 

ntl(i,:)  =  filtf ilt(b,a,ntl(i, : )) ; 
nf 1 (i , : )  =  nt l(i, NzRF+1: end-NzRF) ; 

end 

nfl  =  nf 1/sqrt (var (nfl (:)’)) ; 


7o  zero-pad 

7.  filter 
7o  remove  zeros 

7o  normalize  noise  power 


nt2  =  [zeros (R, NzRF)  n2  zeros (R, NzRF) ] ; 
for  i  =  1:R, 

nt2(i,:)  =  filtf ilt(b,a,nt2(i, : )) ; 
nf2(i,:)  =  nt2 (i , NzRF+1 : end-NzRF) ; 

end 

nf2  =  nf 2/sqrt (var (nf 2 (:)’)) ; 


7.  zero-pad 

7.  filter 
7.  remove  zeros 

7.  normalize  noise  power 


07070707070/070707070707 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7.  Inputs 

07070/0/070/070707070/0/ 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 
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sl_in  =  xfl; 
nl_in  =  nfl; 

Psl_in  =  sum(sl_in . ~2) /length(sl_in) ; 
s2_in  =  xf2; 
n2_in  =  nf2; 

Ps2_in  =  sum(s2_in . ~2) /length(s2_in) ; 


°/„  input  signal  1 

7,  input  noise  1 

7o  power  in  input  signal  1 

*/.  input  signal  2 

7o  input  noise  2 

%  power  in  input  signal  2 


figure(l),  hold  on,  grid 

»/  0  /  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7o  Radiometric  Detection 

07  07  07  07  0  7  0  7  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

w  =  waitbar(O); 

for  i  =  1 : length (SNR) , 

waitbar( i/length (SNR) ,w, [’Radiometer  Progress  -  ’ 
num2str(i) , ’/’ ,num2str (length(SNR) )] ) ; 
nser  =  sqrt (Psl_in/10~ (SNR(i) /10) ) *nl_in;  %  noise  at  required  SNR 
Znr  =  sum (nser . ~2 , 2) ; 

Zsr  =  f lipud(sort (Znr) ) ; 

Tr  =  Zsr(floor (Pfa*R)+l) ; 
for  k  =  1:R, 

Zr(k)  =  sum( (sl_in+nser (k, : ) ) . ~2) ; 

end 

Pdr(i)  =  length(f ind(Zr  >  Tr))/R; 

end 

figure  (1)  ,  plot  (SNR,  Pdr,  ’k-’) ,  xlabel(’SNR  (dB)>),  ylabel( ’P_D’ ) 
close (w) 


°/„  noise  test  statistic 
7,  sort  Zn  (descending) 

7o  find/set  threshold 

%  S+N  test  statistic 

7.  probability  of  detection 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o 

7o  Matched  Filter  Detection 

07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  0 J  07  07  07  07  07  07  0  7  07  07 

/o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

w  =  waitbar(O); 

for  i  =  1 : length (SNR) , 

waitbar( i/length (SNR) ,w, [’Matched  Filter  Progress  - 
num2str(i) ,  ’/’ ,num2str (length(SNR) )] ) ; 
nsem  =  sqrt (Psl_in/10~ (SNR(i) /10) ) *nl_in;  %  noise  at  required  SNR 
for  k  =  1:R, 

Znm(k)  =  sum(sl_in .  *nsem(k,  : ) ) ;  °/0  noise  test  statistic 

end 

Zsm  =  f liplr (sort (Znm) ) ;  %  sort  Zn  (descending) 

Tm  =  Zsm(floor (Pfa*R)+l) ;  %  find/set  threshold 

Zxm  =  sum(sl_in.  ~2) ;  signal  test  statistic 

Zm  =  Zxm  +  Znm;  %  total  test  statistic 


C-ll 


Pdm(i)  =  length(f ind(Zm  >  Tm))/R; 


%  probability  of  detection 


close (w) 

figure (1) ,  plot (SNR, Pdm, ’k— ’) ,  xlabelC’SNR  (dB)’),  ylabel ( ’P_D> ) 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

%  Cross-Correlation  Receiver 

»/  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o 

w  =  waitbar(O);  °/„  create  progress  bar 

for  q  =  1 : length (SNR) , 

nsel  =  sqrt (Psl_in/10~ (SNR(q) /10) ) *nl_in;  %  noise  at  required  SNR 
nse2  =  sqrt (Ps2_in/10~ (SNR(q) /10) ) *n2_in;  %  noise  at  required  SNR 
for  k  =  1:R, 

waitbar (q/length(SNR) ,w, [’CorrRx  1  -  ’ ,num2str(q) , . . . 

’ / ’ ,num2str (length (SNR) ) ,  ’  ;  ' ,num2str (k) , ’ / ’ ,num2str (R)] ) ; 
inpnl  =  nsel(k,:);  %  input  noise  1 

inpsl  =  sl_in  +  nsel(k,:);  °/0  input  signal  1  +  noise  1 


inpn2  =  nse2 (k, : ) ; 
inps2  =  s2_in  +  nse2(k,:); 
outxn  =  xcorr (inpnl , inpn2) 
Zn(k)  =  max(abs (outxn ( : ) ) ) 
outxs  =  xcorr (inpsl , inps2) 
Z(k)  =  max (abs (outxs (:))) ; 


7.  input  noise  2 

7.  input  signal  2  +  noise  2 

l  cross-correlation  (noise) 

7.  max  ’ outxn'  test  statistic 
7.  cross-correlation  (signal+noise) 
7.  max  ’ outxs’  test  statistic 


Zs  =  f liplr (sort (Zn) ) ; 

T  =  Zs (floor (Pfa*R)+l) ; 

Pd(q)  =  length(f ind(Z  >  T))/R; 


7.  sort  test  statistics 
7.  find  threshold 
%  probability  of  detection 


close (w) 

figure (1) ,  plot (SNR, Pd, ’k-o’ ) ,  xlabel(’SNR  (dB)’),  ylabel( JP_D’ ) 


save  chxcorr_wc  SNR  Pdr  Pdm  Pd 


7o  save  required  data  to  workspace 
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%  Brett  D .  Gronholz 
%  EENG  799  —  Sunnner/Fall  2004 
7,  Wideband  Channelized  Receiver 
7.  —  Processing 

0  /  0  /  07  07  0  7  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  0  7  I J  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

clear  all,  close  all,  clc,  format  compact 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  Simulation  Parameters 


07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


Wc  =  le9 ; 

7,  channel  bandwidth 

Nz  =  200; 

7,  one-sided  zero  padding  length  for  ’filtfilt 

fl  =  2 . 5e9 ; 

7,  lower  ChRx  frequency 

fh  =  7 . 5e9 ; 

7,  upper  ChRx  frequency 

o 

II 

o 

T3 

7,  downconvert?  (1  ==  yes,0  ==  no) 

ph  =  0; 

7,  downconversion  starting  phase  (in  degrees) 

meth  =  1 ; 

7,  processing  method  (1=TTM,2=SSM,3=CTM,4=CSM) 

Nfft  =  64; 

7.  FFT  length 

SNR  =  0; 

7.  input  signal  SNR  in  dB 

Nc  =  (fh-f 1) /Wc ; 

7.  number  of  channels 

0/0/ 0/0/ 0/0/ 0/0/ 0/0/ 0/0/ 0,0/ 0/0/ 0/0/ 0/0/ 0/0/ 0/0/ 0/0/0/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  UWB  Signal  Parameters 


07  07  07  07  0  7  0  7  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


fc  =  5e9; 

Tw  =  2/fc; 

Ts  =  2*Tw; 

To  =  Ts/2 ; 
delt  =  0. Ole-9; 
fs  =  1/delt; 

Ns  =  1; 

P  =  1; 
jtr  =  0; 
method  =  ’uni’; 
dly  =  0; 


7.  center  frequency 

7,  pulse  duration 

7o  symbol  duration 

7o  symbol  repetition  interval 

7o  time  resolution 

7o  sample  frequency 

l  number  of  symbols 

7,  signal  power 

7.  jitter  as  percentage  of  Ts 

7.  modulation  method 

7o  first  pulse  delay 


07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 
/o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7o  Generate  Signals 

07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


x  =  uwb (Tw, To, delt, Ns, P, jtr, method, dly) ; 
xs  =  length (x); 

x  =  [zeros (1,240)  x  zeros (1, 240) ] ; 


7.  UWB  signal 
7o  samples/pulse 
7,  add  zeros 
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Px  =  sum(x . ~2) /length(x) ; 
t  =  [0 :delt : (length(x) *delt-delt)] ; 
y  =  sqrt(2*Px)*cos(2*pi*5e9*t) ; 
y  =  [zeros(l,240)  y  zeros(l ,240)] ; 
n  =  randn(l , length (x) ) ; 


7a  power  in  UWB  signal 

7o  time  vector  (length  of  ’x’) 

7,  NB  signal  (tone) 

7,  add  zeros 
7a  generate  noise 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7a  Filter  Input  Signals 

a  /  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  I j  07  07  07  07  07  07 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


N  =  4;  7. 
NzRF  =  200;  7. 
[b,a]  =  butter (N, [f l/(f s/2)  fh/(fs/2)] ) ; % 


filter  order 

one-sided  zero  padding  length 
filter  coeffs 


xf  =  [zeros (1, NzRF)  x  zeros(l,NzRF)] ;  % 

xf  =  real (f  iltf  ilt  (b, a,xf ) ) ;  7, 

xf  =  xf  (NzRF+1 :  end-NzRF)  ;  7. 

Pxf  =  sum  (xf  .  ~  2) /length  (xf) ;  7, 

yf  =  [zeros (1, NzRF)  y  zeros(l,NzRF)]  ;  7o 
yf  =  real (f iltf ilt (b, a,yf ) ) ;  7o 

yf  =  yf  (NzRF+1:  end-NzRF)  ;  7. 

yf  =  sqrt(Pxf/var(yf  ))*yf  ;  7« 

nf  =  [zeros (1, NzRF)  n  zeros(l,NzRF)]  ;  7o 
nf  =  filtf ilt(b,a,nf ) ;  7o 

nf  =  nf  (NzRF+1:  end-NzRF)  ;  7. 

nf  =  nf /sqrt  (var  (nf ) ) ;  7o 

nf  =  sqrt  (Pxf /10~  (SNR/10) )  *nf ;  7« 


07  07  07  07  0  7  0  7  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  0  7  0  7  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%  Channelized  Receiver 


07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  0  7  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


inp  =  xf ; 

out  =  chrx(inp,fs,Wc,fl,fh,Nz,dc,ph) ; 


07  07  07  07  0  7  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  I j 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%  Channelized  Receiver  Processing 

07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  0 j  0  7  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


if  meth  ==  1, 

outx  =  fftshift(ifft(out,Nfft, 1) , 1) ; 
elseif  meth  ==  2, 

outx  =  fftshift(fft(out,Nfft,2) ,2) ; 


zero-pad 

filter 

remove  zeros 

power  in  filtered  UWB  signal 

zero-pad 

filter 

remove  zeros 
scale 

zero-pad 

filter 

remove  zeros 
normalize  noise  power 
noise  at  required  SNR 


7.  input  sig 
7,  channelized  receiver 

7,  Temp-Temp  Matrix 
7,  Spec-Spec  Matrix 
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elseif  meth  ==  3, 

outx  =  ctm(out ,Nf ft , 1 , 0) ;  %  Cross-Temporal  Matrix 

elseif  meth  ==  4, 

outx  =  csm(out ,Nf ft , 1 , 0) ;  %  Cross-Spectral  Matrix 

else 

error ( ' Invalid  value  (meth)’) 

end 

figure (1) ,pcolor (abs (outx) ) , shading ( ' interp' ) , colormap( 1-gray) 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


°/„  Sum  of  Real/Imaginary  Parts  of  Processed  Matrix 

»/  07  07  07  07  07  07  07  07  07  07  07  07  0 j  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0 j  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07  I j  0  7  07  07 

/o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0 

figure (2) 

subplot (2 , 2 , 1) , plot (real (sum(outx, 1) ), ’k' ) ,title( 'Real  -  Sum  of  Rows') 
axis([0  length (outx(l ,:) )  ... 

min(real(sum(outx, l)))-le-4  max(real(sum(outx, l)))+le-4] ) 
subplot (2 , 2 , 2) , plot (imag( sum (outx, 1) ) , 'k' ) ,title( ' Imaginary  -  Sum  of  Rows') 
axis([0  length (outx(l ,:) )  ... 

min (imag( sum (outx, l)))-le-4  max (imag( sum (outx, l)))+le-4] ) 
subplot (2 , 2 , 3) ,plot (real (sum(outx, 2) ) , 'k' ) , title( 'Real  -  Sum  of  Columns’) 
axis([0  Nfft  min(real (sum(outx , 2) ) )-le-4  max(real(sum(outx,2)))+le-4] ) 
subplot (2 , 2 ,4) , plot (imag( sum (outx, 2) ), 'k' ) ,title( ' Imaginary  -  Sum  of  Columns') 
axis([0  Nfft  min (imag( sum (outx , 2) ) )-le-4  max(imag(sum(outx,2)))+le-4] ) 


07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7o  Plot  Channelized  Receiver  Outputs 

07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o 


figured) 
for  i  =  l:Nc, 

subplot (Nc , 1 , i) ,plot (out (i , : ) , 'k' ) 

axis([0  length(out)  min(out ( : ) ) -le-1  max(out( : ) )  +  le-l] ) 
ylabel(['Ch  ' ,num2str (i)] ) 

end 
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%  Brett  D .  Gronholz 
%  EENG  799  —  Sunnner/Fall  2004 
%  UWB  Detection  Probability  (Pd) 

%  —  Varying  Points  in  (I)FFT  and  SNR 

y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o 

clear  all,  close  all,  clc,  format  compact 


0/  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7o  Simulation  Parameters 

07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


Wc  = 

[le9] ; 

7.  channel  bandwidth  to  simulate 

Nz  = 

[200] ; 

7.  zero-padding  length  for  ’Wc’ 

fl  = 

2 . 5e9 ; 

7.  lower  ChRx  frequency 

fh  = 

7 . 5e9 ; 

7.  upper  ChRx  frequency 

dc  = 

0; 

7.  downconvert?  (l=yes,0=no) 

ph  = 

0; 

7.  DC  mixer  starting  phase 

meth 

=  1; 

7,  processing  method  (1=TTM,2=SSM,3=CTM,4=CSM) 

Nfft  =  [4,8,16,32,64,128,256];  7.  FFT  lengths  to  simulate  (TTM/CTM) 

%  Nfft  =  [512,1024,  2048];  7,  FFT  lengths  to  simulate  (SSM/CSM) 

SNR  =  [-10:1:5];  %  SNR  to  simulate 

Pfa  =  10~-2;  7.  probability  of  false  alarm 

R  =  10/Pfa;  7.  number  of  realizations 


07  07  07  07  0  7  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  UWB  Signal  Parameters 

07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


fc  =  5e9; 

Tw  =  2/fc; 

Ts  =  2*Tw; 

To  =  Ts/2 ; 
delt  =  0. Ole-9; 
fs  =  1/delt; 

Ns  =  1; 

P  =  1; 
jtr  =  0; 
method  =  ’uni’; 
dly  =  0; 


7.  center  frequency 

7o  pulse  duration 

7.  symbol  duration 

7o  symbol  repetition  interval 

7o  time  resolution 

7o  sample  frequency 

7,  number  of  symbols 

7«  signal  power 

7.  jitter  as  percentage  of  Ts 

7o  modulation  method 

7.  first  pulse  delay 


07  07  07  07  07  07  07  07  07  07  07  07  I J  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7.  Generate  Signals 

07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


x  =  uwb (Tw, To, delt, Ns, P, jtr, method, dly) ; 


7.  UWB  signal 
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x  =  [zeros (1 ,240)  x  zeros (1, 240) ] ; 
Px  =  sum(x . ~2) /length(x) ; 
t  =  [0 :delt : (length(x)*delt-delt)] ; 
y  =  sqrt(2*Px)*cos(2*pi*5e9*t) ; 
n  =  randn(R, length (x) ) ; 


l  add  zeros 
%  power  in  UWB  signal 
%  time  vector  (length  of  Jx’) 
7,  NB  signal  (tone) 

7a  ’  R;  noise  realizations 


0/  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%  Filter  Input  Signals  (Hrf) 

07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  0 J  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

N  =  4;  7,  order  of  Hrf  BPF 

NzRF  =  200;  °/0  one-sided  zero  padding  length 

[b,a]  =  butter (N,  [fl/(f s/2)  fh/(fs/2)] ) ; %  Hrf  filter  coeffs 


xf  =  [zeros (1, NzRF)  x  zeros(l,NzRF)] ; 
xf  =  real (f iltf ilt (b, a,xf ) ) ; 
xf  =  xf (NzRF+1 : end-NzRF) ; 

Pxf  =  sum(xf . ~2) /length (xf) ; 

yf  =  [zeros (1, NzRF)  y  zeros(l,NzRF)] ; 
yf  =  real (f iltf ilt (b, a,yf ) ) ; 
yf  =  yf (NzRF+1: end-NzRF) ; 
yf  =  sqrt(Pxf/var(yf ))*yf ; 

nt  =  [zeros (R, NzRF)  n  zeros(R,NzRF)] ; 
for  i  =  1:R, 

nt(i,:)  =  f iltf ilt (b , a,nt (i , : ) ) ; 
nf(i,:)  =  nt(i, NzRF+1 : end-NzRF) ; 

end 

nf  =  nf /sqrt (var(nf (:)’)) ; 


7o  zero-pad 
7,  filter 
7,  remove  zeros 

7.  power  in  filtered  UWB  signal 

7o  zero-pad 
7o  filter 
7o  remove  zeros 

7o  scale  (power  equal  to  ’Pxf’) 

7o  zero-pad 

7o  filter 
7o  remove  zeros 

7o  normalize  noise  power 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 
/o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  ChRx  Inputs 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o 


s_in  =  xf ;  %  input  signal 

n_in  =  nf ;  %  input  noise 

Ps_in  =  sum(s_in . ~2) /length(s_in) ; 


7,  power  in  input  signal 


figure(l),  hold  on,  grid 
for  c  =  1 : length (Nf f t) , 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 
/o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%  Channelized  Receiver 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0 


7,  Nfft  loop 
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w  =  waitbar(O);  %  create  progress  bar 

for  q  =  1 : length (SNR) ,  7,  SNR  loop 

nse  =  sqrt  (Ps_in/10~  (SNR(q) /10) ) *n_in;  7„  noise  at  required  SNR 
for  k  =  1:R,  %  realizations  loop 

waitbar (q/length(SNR) ,w, [’ChRx  Progress  -  ’ ,num2str (q) , ’ / ’ , . . . 

num2str (length(SNR) ) , ’ ;  ’ ,num2str(k) ,  ’/’ ,num2str(R)] ) ; 
inpn  =  nse(k,:);  °/„  input  N 

inps  =  s_in  +  nse(k,:);  %  input  S+N 

outn  =  chrx(inpn,fs,Wc,fl,fh,Nz,dc,ph) ;  %  output  N 
outs  =  chrx(inps,fs,Wc,fl,fh,Nz,dc,ph) ;  °/„  output  S+N 
if  meth  ==  1 , 

outxn  =  ifft(outn,Nfft(c) , 1) ;  7,  Temp-Temp  Matrix,  N 

outxs  =  ifft(outs,Nfft(c) ,  1)  ;  7,  Temp-Temp  Matrix,  S+N 

elseif  meth  ==  2, 

outxn  =  fft(outn,Nfft(c)  ,2)  ;  7„  Spec-Spec  Matrix,  N 

outxs  =  fft(outs,Nfft(c)  ,2)  ;  7  Spec-Spec  Matrix,  S+N 

elseif  meth  ==  3, 

outxn  =  ctm(outn,Nfft(c) ,  1,0)  ;  7  Cross-Temporal  Matrix,  N 

outxs  =  ctm(outs,Nfft(c) ,  1,0)  ;  7  Cross-Temporal  Matrix,  S+N 

elseif  meth  ==  4, 

outxn  =  csm(outn,Nfft(c) ,  1,0)  ;  7  Cross-Spectral  Matrix,  N 

outxs  =  csm(outs,Nfft(c) ,  1,0)  ;  7  Cross-Spectral  Matrix,  S+N 

else 

error (’ Invalid  value  (meth)’) 

end 

Zn(k)  =  max(abs (outxn( : ) ) ) ;  /  max  ’outxn’  test  statistic 

Z(k)  =  max (abs (outxs (:))) ;  7  max  ’outxs’  test  statistic 

end 

Zs  =  f liplr (sort (Zn) ) ;  /  sort  ’outxn’  test  statistics 

T  =  Zs(f loor (Pf a*R)+l)  ;  7,  find/set  threshold 

Pd(c,q)  =  length (f ind(Z  >  T))/R;  70  probability  of  detection 

end 

close (w) 
if  c  ==  1 , 

figure (1) ,  plot (SNR,Pd(c , : ) , ’k-o ’ ) ,  xlabel(’SNR  (dB)’),  ylabel(’P_D’) 
elseif  c  ==  2, 

figure (1) ,  plot (SNR,Pd(c , : ) , ’k-s ’ ) ,  xlabel(’SNR  (dB)’),  ylabel ( ’P_D ’ ) 
elseif  c  ==  3, 

figure (1) ,  plot (SNR,Pd(c , : ) , ’k-~ ’ ) ,  xlabel(’SNR  (dB)’),  ylabel ( ’P_D ’ ) 
elseif  c  ==  4, 

figure (1) ,  plot (SNR,Pd(c , : ) , ’k-d’ ) ,  xlabel(’SNR  (dB)’),  ylabel(’P_D’) 
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elseif  c  ==  5, 

figure (1) ,  plot (SNR,Pd(c , : ) , ’k-* 1 ) ,  xlabel(’SNR  (dB)’),  ylabel ( >P_D 1 ) 
elseif  c  ==  6, 

figure (1) ,  plot (SNR,Pd(c , : ) , ,k-+) ) ,  xlabel(’SNR  (dB)'),  ylabel ( >P_D 1 ) 
elseif  c  ==  7, 

figure (1) ,  plot (SNR,Pd(c , : ) , ’k-x’ ) ,  xlabel(’SNR  (dB)'),  ylabel ( ’P_D' ) 

end 

end 

save  ttml_nfft  SNR  Pd  7„  save  required  data  to  workspace 
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%  Brett  D .  Gronholz 
%  EENG  799  —  Sunnner/Fall  2004 
%  UWB  Detection  Probability  (Pd) 

%  —  Varying  Channel  Bandwidth  and  SNR 

y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o 

clear  all,  close  all,  clc,  format  compact 


0/  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7o  Simulation  Parameters 

07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

Wc  =  [Ie9,0.5e9,0.25e9,0. Ie9] ;  %  channel  bandwidths  to  simulate 

Nz  =  [200,400,800,2000];  °/0  zero-padding  lengths  for  each  ’Wc’ 


fl  =  2 . 5e9 ; 
fh  =  7 . 5e9 ; 
dc  =  0; 
ph  =  0; 
meth  =  1 ; 

Nfft  =  64; 

SNR  =  [-10:1:5] ; 
Pfa  =  10~-2 ; 

R  =  10/Pfa; 


%  lower  Hrf/ChRx  frequency 
%  upper  Hrf/ChRx  frequency 
7.  downconvert?  (l=yes,0=no) 

7.  DC  phases  to  simulate 

7,  processing  method  (1=TTM,2=SSM,3=CTM,4=CSM) 
7.  (I)FFT  length 
7.  SNRs  to  simulate 
7.  probability  of  false  alarm 
7.  number  of  realizations 


07  07  07  07  07  07  07  07  07  07  07  07  07  0 J  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  UWB  Signal  Parameters 


07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


fc  =  5e9; 

Tw  =  2/fc; 

Ts  =  2*Tw; 

To  =  Ts/2 ; 
delt  =  0. Ole-9; 
fs  =  1/delt; 

Ns  =  1; 

P  =  1; 
jtr  =  0; 
method  =  ’uni’; 
dly  =  0; 


7.  center  frequency 
7o  pulse  duration 
7o  symbol  duration 
7.  symbol  repetition  interval 
7o  time  resolution 
7o  sample  frequency 
7,  number  of  symbols 
I  signal  power 

7o  jitter  as  percentage  of  Ts 
7o  modulation  method 
7«  first  pulse  delay 


07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  Generate  Signals 

07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  0 J  0  7  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


x  =  uwb (Tw, To, delt, Ns, P, jtr, method, dly) ; 
x  =  [zeros (1,240)  x  zeros (1, 240) ] ; 


7,  UWB  signal 
7,  add  zeros 
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Px  =  sum(x . ~2) /length(x) ; 
t  =  [0 :delt : (length(x)*delt-delt)] ; 
y  =  sqrt(2*Px)*cos(2*pi*5e9*t) ; 
n  =  randn(R, length (x) ) ; 


%  power  in  UWB  signal 

7o  time  vector  (length  of  ’x’) 

7o  NB  signal  (tone) 

7a  ’IT  noise  realizations 


0/  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7o  Filter  Input  Signals 

07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


N  =  4;  7,  order  of  Hrf  BPF 

NzRF  =  200;  °/0  one-sided  zero  padding  length 

[b,a]  =  butter(N,  [f l/(f s/2)  fh/(fs/2)] ) Hrf  filter  coeffs 


xf  =  [zeros (1, NzRF)  x  zeros(l,NzRF)] ; 
xf  =  real (f iltf ilt (b, a,xf ) ) ; 
xf  =  xf (NzRF+1 : end-NzRF) ; 

Pxf  =  sum(xf . ~2) /length (xf) ; 

yf  =  [zeros (1, NzRF)  y  zeros(l,NzRF)] ; 
yf  =  real (f iltf ilt (b, a,yf ) ) ; 
yf  =  yf (NzRF+1: end-NzRF) ; 
yf  =  sqrt(Pxf/var(yf))*yf ; 


7o  zero-pad 
7,  filter 
7.  remove  zeros 

7.  power  in  filtered  UWB  signal 

7.  zero-pad 
7o  filter 
7o  remove  zeros 
7o  scale 


nt  =  [zeros (R, NzRF)  n  zeros (R, NzRF)] ; 
for  i  =  1:R, 

nt(i,:)  =  f iltf ilt (b , a,nt (i , : ) ) ; 
nf(i,:)  =  nt(i, NzRF+1 : end-NzRF) ; 

end 

nf  =  nf /sqrt (var(nf (:)’)) ; 


7o  zero-pad 

7o  filter 
7o  remove  zeros 

7.  normalize  noise  power 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  ChRx  Inputs 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


s_in  =  xf ;  %  input  signal 

n_in  =  nf ;  %  input  noise 

Ps_in  =  sum(s_in . ~2) /length(s_in) ; 


7,  power  in  input  signal 


figure(l),  hold  on,  grid 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%  Radiometric  Detection 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


w  =  waitbar(O);  %  create  progress  bar 

for  i  =  1 :  length  (SNR) ,  “/,  SNR  loop 
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waitbar( i/length (SNR) ,w, [’ Radiometer  Progress  -  ’ ,  ... 

num2str(i) ,  ’/’ ,num2str (length (SNR) )] ) ; 
nser  =  sqrt(Ps_in/10~(SNR(i)/10))*n_in;  “/„  noise  at  required  SNR 
Znr  =  sum(nser . ~2 , 2) ;  %  noise  test  statistic 

Zsr  =  f lipud(sort (Znr) ) ;  %  sort  Zn  (descending) 

Tr  =  Zsr(floor (Pfa*R)+l) ;  %  find/set  threshold 

for  k  =  1:R, 


Zr(k)  =  sum( (s_in+nser (k, : ) ) . ~2) ;  %  S+N  test  statistics 


end 


Pdr(i)  =  length(f ind(Zr  >  Tr))/R;  %  probability  of  detection 

end 

close (w) 

figure  (1)  ,  plot  (SNR,  Pdr,  ’k-’ ) ,  xlabel(’SNR  (dB)’),  ylabel( ’P_D’ ) 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%  Matched  Filter  Detection 


0  /  07  07  07  0  7  I J  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07 

/o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


w  =  waitbar(O);  %  create  progress  bar 

for  i  =  1 : length (SNR) ,  %  SNR  loop 

waitbar( i/length (SNR) ,w, [’Matched  Filter  Progress  - 
num2str(i) ,  ’/’ ,num2str (length (SNR) )] ) ; 
nsem  =  sqrt(Ps_in/10~(SNR(i)/10))*n_in;  “/„  noise  at  required  SNR 
for  k  =  1:R, 

Znm(k)  =  sum(s_in . *nsem(k, : ) ) ;  %  noise  test  statistic 


end 

Zsm  =  f liplr (sort (Znm) ) ; 

Tm  =  Zsm(floor (Pfa*R)+l) ; 

Zxm  =  sum(s_in . ~2) ; 

Zm  =  Zxm  +  Znm; 

Pdm(i)  =  length(f ind(Zm  >  Tm))/R; 

end 

close (w) 

figure  (1)  ,  plot  (SNR,  Pdm,  ’k— ’) ,  xlabel(’SNR 


%  sort  Znm  (descending  value) 
*/.  find  threshold  based  on  Pfa 
%  signal  test  statistic 
7o  total  test  statistic 
7.  probability  of  detection 


(dB)’),  ylabel ( ’P_D ’ ) 


for  c  =  1 :  length(Wc) ,  °/0  ChBW  loop 

07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

7o  Channelized  Receiver 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


w  =  waitbar(O);  °/0  create  progress  bar 

for  q  =  1 : length (SNR) ,  %  SNR  loop 

nse  =  sqrt (Ps_in/10~ (SNR(q) /10) ) *n_in;  %  noise  at  required  SNR 
for  k  =  1:R,  °/„  realizations  loop 
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waitbar (q/length(SNR) ,w, [’ChRx  Progress  -  ’ ,num2str (q) , ’ / ’ , . . . 

num2str (length(SNR) ) , ’ ;  ’ ,num2str(k) ,  ’/’ ,num2str(R)] ) ; 
inpn  =  nse(k,:);  °/„  input  N 

inps  =  s_in  +  nse(k,:);  °/„  input  S+N 

outn  =  chrx(inpn,fs,Wc(c)  ,fl,fh,Nz(c)  ,dc,ph) ;  °/0  output  N 
outs  =  chrx(inps,fs,Wc(c)  ,fl,fh,Nz(c)  ,dc,ph) ;  °/0  output  S+N 
if  meth  ==  1, 

outxn  =  ifft(outn,Nfft,  1)  ;  °/0  Temp-Temp  Matrix,  N 

outxs  =  ifft(outs,Nfft,  1)  ;  °/0  Temp-Temp  Matrix,  S+N 

elseif  meth  ==  2, 

outxn  =  fft(outn,Nfft,2)  ;  °/0  Spec-Spec  Matrix,  N 

outxs  =  fft(outs,Nfft,2) ;  °/0  Spec-Spec  Matrix,  S+N 

elseif  meth  ==  3, 

outxn  =  ctm(outn,Nf ft ,  1 , 0) ;  °/0  Cross-Temporal  Matrix,  N 
outxs  =  ctm(outs ,Nf ft ,  1 , 0) ;  °/0  Cross-Temporal  Matrix,  S+N 
elseif  meth  ==  4, 

outxn  =  csm(outn,Nf ft ,  1 , 0) ;  °/0  Cross-Spectral  Matrix,  N 
outxs  =  csm(outs ,Nf ft ,  1 , 0) ;  °/0  Cross-Spectral  Matrix,  S+N 

else 

error (’ Invalid  value  (meth)’) 

end 

Zn(k)  =  max(abs (outxn( : ) ) )  ;  °/0  max  ’outxn’  test  statistic 

Z(k)  =  max (abs (outxs (:)))  ;  °/0  max  ’outxs’  test  statistic 

end 

Zs  =  f liplr (sort (Zn) )  ;  °/0  sort  ’outxn’  test  statistics 

T  =  Zs(f loor (Pf a*R)+l)  ;  °/0  find/set  threshold 

Pd(c,q)  =  length(f ind(Z  >  T))/R;  °/0  probability  of  detection 

end 

close (w) 
if  c  ==  1 , 

figure (1) ,  plot (SNR,Pd(c , : ) , ’k-o ’ ) ,  xlabel(’SNR  (dB)’),  ylabel ( ’P_D ’ ) 
elseif  c  ==  2, 

figure (1) ,  plot (SNR,Pd(c , : ) , ’k-s ’ ) ,  xlabel(’SNR  (dB)’),  ylabel ( ’P_D ’ ) 
elseif  c  ==  3, 

figure (1) ,  plot (SNR,Pd(c , : ) , ’k-“ ’ ) ,  xlabel(’SNR  (dB)’),  ylabel( ’P_D’ ) 
elseif  c  ==  4, 

figure (1) ,  plot (SNR,Pd(c , : ) , ’k-d’ ) ,  xlabel(’SNR  (dB)’),  ylabel ( ’P_D ’ ) 

end 

end 

save  ttm64_wc  SNR  Pdr  Pdm  Pd  °/0  save  required  data  to  workspace 
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%  Brett  D .  Gronholz 
%  EENG  799  —  Sunnner/Fall  2004 
%  UWB  Detection  Probability  (Pd) 

%  —  Varying  Downconversion  Mixer  Phase  and  Channel  Bandwidth 

%  (fixed  SNR) 


y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


clear  all,  close  all,  clc,  format  compact 


»/  07  07  07  07  0  7  07  07  07  07  07  07  07  0 J  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  Simulation  Parameters 


07  07  07  07  07  0  7  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o 


Wc  = 
Nz  = 
fl  = 
fh  = 
dc  = 
ph  = 
meth 
Nfft 
SNR 
Pfa 
R  = 


[le9]  ; 

7.  channel  bandwidths  to  simulate 

[200]  ; 

7.  zero-padding  lengths  for  each  ’Wc’  (one-sic 

2 . 5e9 ; 

7.  lower  ChRx  frequency 

7 . 5e9 ; 

7.  upper  ChRx  frequency 

l; 

7.  downconvert?  (l=yes,0=no) 

[0:1:360] ; 

7.  phases  to  simulate 

=  l; 

7,  processing  method  (1=TTM,2=SSM,3=CTM,4=CSM) 

=  64; 

7.  FFT  length 

=  0; 

7.  SNR  to  simulate 

=  10~-2 ; 

7.  probability  of  false  alarm 

10/Pfa; 

7.  number  of  realizations 

07  07  07  07  0  7  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07 

/o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  UWB  Signal  Parameters 

07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


fc  =  5e9; 

Tw  =  2/fc; 

Ts  =  2*Tw; 

To  =  Ts/2 ; 
delt  =  0. Ole-9; 
fs  =  1/delt; 

Ns  =  1; 

P  =  1; 
jtr  =  0; 
method  =  ’uni’; 
dly  =  0; 


7.  center  frequency 

7o  pulse  duration 

7.  symbol  duration 

7o  symbol  repetition  interval 

7o  time  resolution 

7o  sample  frequency 

7,  number  of  symbols 

7«  signal  power 

7.  jitter  as  percentage  of  Ts 

7o  modulation  method 

7.  first  pulse  delay 


07  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  0  7  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o 


7.  Generate  Signals 

07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


x  =  uwb (Tw, To, delt, Ns, P, jtr, method, dly) ; 


7,  UWB  signal 
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x  =  [zeros (1 ,240)  x  zeros (1, 240) ] ; 
Px  =  sum(x . ~2) /length(x) ; 
t  =  [0 :delt : (length(x)*delt-delt)] ; 
y  =  sqrt(2*Px)*cos(2*pi*5e9*t) ; 
n  =  randn(R, length (x) ) ; 


l  add  zeros 
%  power  in  UWB  signal 
%  time  vector  (length  of  Jx’) 
7,  NB  signal  (tone) 

7a  ’  R;  noise  realizations 


0/  07  07  07  07  07  07  07  07  07  07  07  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%  Filter  Input  Signals  (Hrf) 

07  07  07  07  07  07  07  07  07  07  07  07  0  7  0  7  07  07  07  07  07  07  07  07  07  07  07  07  07  0 J  07  07  07  07 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

N  =  4;  7,  order  of  Hrf  BPF 

NzRF  =  200;  °/0  one-sided  zero  padding  length  for 

[b,a]  =  butter (N,  [fl/(f s/2)  fh/(fs/2)] ) ; %  filter  coeffs 


’filtfilt’ 


xf  =  [zeros (1, NzRF)  x  zeros(l,NzRF)] ; 
xf  =  real (filtfilt (b, a, xf) ) ; 
xf  =  xf (NzRF+1 : end-NzRF) ; 

Pxf  =  sum(xf . ~2) /length (xf) ; 

yf  =  [zeros (1, NzRF)  y  zeros(l,NzRF)]  ; 
yf  =  real (filtfilt (b, a, yf) ) ; 
yf  =  yf (NzRF+1: end-NzRF) ; 
yf  =  sqrt(Pxf/var(yf ))*yf ; 


7o  zero-pad 
7,  filter 
7,  remove  zeros 

7.  power  in  filtered  UWB  signal 

7o  zero-pad 
7o  filter 
7o  remove  zeros 
7o  scale 


nt  =  [zeros (R, NzRF)  n  zeros(R,NzRF)] ; 
for  i  =  1:R, 

nt(i,:)  =  f iltf ilt (b , a,nt (i , : ) ) ; 
nf(i,:)  =  nt(i, NzRF+1 : end-NzRF) ; 

end 

nf  =  nf /sqrt (var(nf (:)’)) ; 


7o  zero-pad 

7o  filter 
7o  remove  zeros 

7o  normalize  noise  power 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 
/o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7o  ChRx  Inputs 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 

/o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o 


s_in  =  xf ;  %  input  signal 

n_in  =  nf ;  %  input  noise 

Ps_in  =  sum(s_in . ~2) /length(s_in) ; 


7,  power  in  input  signal 


figure(l),  hold  on,  grid 
for  c  =  1 : length (Wc) , 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 
/o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%  Channelized  Receiver 


07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07  07 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


7,  ChBW  loop 
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w  =  waitbar(O);  7„  create  progress  bar 

nse  =  sqrt  (Ps_in/10~  (SNR/10) )  *n_in;  °/„  noise  at  required  SNR 

for  q  =  1 : length(ph) ,  %  phase  loop 

for  k  =  1:R,  %  realizations  loop 

waitbar (q/length(ph) ,w, [’ChRx  1  -  ’ ,num2str(q) . 

num2str (length(ph) ) , ’ ;  ’ ,num2str(k) ,  ’/’ ,num2str(R)] ) ; 
inpn  =  nse(k,:);  °/„  input  N 

inps  =  s_in  +  nse(k,:);  °/„  input  S+N 

outn  =  chrx(inpn,fs,Wc(c) ,fl,fh,Nz(c) ,dc,ph(q)) ;  %  output  N 

outs  =  chrx(inps,fs,Wc(c) ,fl,fh,Nz(c) ,dc,ph(q)) ;  %  output  S+N 

if  meth  ==  1 , 

outxn  =  ifft(outn,Nfft, 1) ;  °/0  Temp-Temp  Matrix,  N 

outxs  =  ifft(outs,Nfft,  1)  ;  °/0  Temp-Temp  Matrix,  S+N 

elseif  meth  ==  2, 

outxn  =  fft(outn,Nfft,2) ;  °/0  Spec-Spec  Matrix,  N 

outxs  =  fft(outs,Nfft,2) ;  °/0  Spec-Spec  Matrix,  S+N 

elseif  meth  ==  3, 

outxn  =  ctm(outn,Nf ft ,  1 , 0) ;  °/0  Cross-Temporal  Matrix,  N 
outxs  =  ctm(outs ,Nf ft ,  1 , 0) ;  °/0  Cross-Temporal  Matrix,  S+N 
elseif  meth  ==  4, 

outxn  =  csm(outn,Nf ft ,  1 , 0) ;  °/0  Cross-Spectral  Matrix,  N 
outxs  =  csm(outs ,Nf ft ,  1 , 0) ;  °/0  Cross-Spectral  Matrix,  S+N 

else 

error (’ Invalid  value  (meth)’) 

end 

Zn(k)  =  max(abs (outxn( : ) ) )  ;  °/0  max  ’outxn’  test  statistic 

Z(k)  =  max (abs (outxs (:)))  ;  °/0  max  ’outxs’  test  statistic 

end 

Zs  =  f liplr (sort (Zn) )  ;  °/0  sort  ’outxn’  test  statistics 

T  =  Zs(f loor (Pf a*R)+l)  ;  °/0  find/set  threshold 

Pd(c,q)  =  length(f ind(Z  >  T))/R;  °/0  probability  of  detection 

end 

close (w) 
if  c  ==  1 , 

figure(l),  plot (ph,Pd(c , : ) , ’k-o ’ ) ,  xlabel( ’Phase  (deg)’),  ylabel ( ’P_D ’ ) 
elseif  c  ==  2, 

figure(l),  plot (ph,Pd(c, : ) , ’k-s ’ ) ,  xlabel( ’Phase  (deg)’),  ylabel ( ’P_D ’ ) 
elseif  c  ==  3, 

figure(l),  plot (ph,Pd(c, : ) , ’k-~ ’ ) ,  xlabel( ’Phase  (deg)’),  ylabel ( ’P_D ’ ) 
elseif  c  ==  4, 

figure(l),  plot (ph,Pd(c, : ) , ’k-d’ ) ,  xlabel( ’Phase  (deg)’),  ylabel ( ’P_D ’ ) 

end 
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end 


save  ttml_ph  ph  Pd 


%  save  required  data  to  workspace 


C-27 


function  sig  =  uwb(Tw, To, delt, Ns, P,jtr, method, dly) ; 

7oUWB  Ultra-Wideband  Signal  Generator 
%  SIG  =  UWB(Tw, To, delt, Ns, P,jtr, method, dly) 

1 

7.  INPUTS 

%  Tw  -  pulse  duration 

l  To  -  symbol  repetition  interval 

%  delt  -  time  resolution 

%  Ns  -  number  of  symbols 

%  P  -  signal  power 

7.  jtr  -  jitter  as  percentage  of  Ts  =  2*Tw 

%  method  -  'uni',  ’ppm',  'pam' ,  or  'bppm' 

%  dly  -  first  pulse  delay 

% 

7a  OUTPUT 

7a  sig  -  UWB  output  signal 

7. 

1 

7a  Brett  D.  Gronholz 

7a  EENG  799  —  Summer/Fall  2004 

7. 

if  nargin  ~=  8, 

error ('Not  enough  input  arguments!') 

end 

if  delt  <=  0, 

error('Time  resolution  must  be  greater  than  0!’); 

end 

if  Ns  <  1, 

error ( 'Number  of  symbols  must  be  greater  than  zero!’) 

end 

if  P  <=  0, 

error ('Power  must  be  greater  than  zero!') 

end 

if  jtr  >  1  I  jtr  <  0, 

error (' Jitter  percentage  must  be  between  0  and  1.') 

end 

if  dly  <  0, 

error('First  pulse  delay  time  must  be  greater  than  or  equal  to  0.') 

end 
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rand(' state' , sum(100*clock) ) 


%  Variables 

Tm  =  0.4*Tw; 

t  =  0 : delt : (2*Tw-delt) ; 

Nc  =  length(t) ; 

Nr  =  length(0 : delt : To-delt) ; 


°/„  pulse  width  parameter 
/  time  vector  for  UWB  pulse  generation 
/  number  of  samples/symbol 
%  number  of  samples/repetition  interval 


%  Generate  UWB  pulse  (2nd  derivative  of  Gaussian  pulse) 

w  =  (l-4*pi* ( (t-1 . 2*Tm) /Tm) . "2) . *exp(-2*pi* ( (t — 1 . 2*Tm) /Tm) . "2) ;  °/„  UWB  pulse 
Pw  =  (1/To) *sum(w . ~2) *delt ;  %  power  in  w 

s  =  sqrt(P/Pw)*w;  %  received  UWB  waveform 


if  strcmpi (method, 'uni ’) , 
sig  =  zeros (1 ,Ns*Nr) ; 
r  =  round(2*(rand(l,Ns)-0.5)*jtr*Nc/2) ; 
r (1)  =  0; 
for  i  =  l:Ns, 

sig(Nr*(i-l)+l+r(i) :Nr*(i-l)+Nc/2+r(i) )  =  s(l:Nc/2); 

end 


elseif  strcmpi (method, 'ppm' ) , 
ppml  =  s ; 

ppmO  =  fliplr(ppml) ; 
sig  =  zeros ( 1, Ns*Nr) ; 
r  =  round(2*(rand(l,Ns)-0.5)*jtr*Nc) ; 
r (1)  =  0; 

bits  =  randint (1 ,Ns)  ; 
for  i  =  l:Ns, 

if  bits(i)  ==  0, 

sig(Nr*(i-l)+l+r(i) :Nr*(i-l)+Nc+r(i))  =  ppmO; 

else 

sig(Nr*(i-l)+l+r(i) :Nr*(i-l)+Nc+r(i))  =  ppml; 

end 

end 

elseif  strcmpi (method, 'pam' ) , 
paml  =  sqrt (0 . 5) *s (1 : Nc/2) ; 
pamO  =  -paml ; 
sig  =  zeros (1 ,Ns*Nr) ; 

r  =  round (2* (rand(l , Ns)-0 . 5) *jtr*Nc/2) ; 
r (1)  =  0; 
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bits  =  randint (1 ,Ns) ; 
for  i  =  l:Ns, 

if  bits(i)  ==  0, 

sig(Nr*(i-l)+l+r(i) :Nr*(i-l)+Nc/2+r(i))  = 

else 

sig(Nr*(i-l)+l+r(i) :Nr*(i-l)+Nc/2+r(i))  = 

end 

end 

elseif  strcmpi (method, ’bppm' ) , 
bppmlO  =  s; 
bppmOl  =  -s; 

bppmOO  =  fliplr (bppmOl) ; 

bppmll  =  fliplr (bppmlO) ; 

sig  =  zeros (1 ,Ns*Nr) ; 

r  =  round(2*(rand(l,Ns)-0.5)*jtr*Nc) ; 

r (1)  =  0; 

bits  =  randint ( 1, 2*Ns) ; 
for  i  =  l:Ns, 

switch  bi2de ( [bits (2*i-l)  bits(2*i)]) 
case  [0] 

sig(Nr*(i-l)+l+r(i) :Nr*(i-l)+Nc+r(i)) 
case  [1] 

sig(Nr* (i-l)+l+r (i) :Nr*(i-l)+Nc+r(i)) 
case  [2] 

sig(Nr* (i-l)+l+r (i) :Nr*(i-l)+Nc+r(i)) 
case  [3] 

sig(Nr*(i-l)+l+r(i) :Nr*(i-l)+Nc+r(i)) 

end 

end 


else 

error ( 'Modulation  type  error.'); 

end 

sig  =  [zeros (1 , round (dly/delt) )  sig]; 


pamO ; 
paml ; 


=  bppmOO; 
=  bppmO 1 ; 
=  bppmlO; 
=  bppmll; 
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function  out  =  chrx(inp,fs,Wc,fl,fh,Nz,dc,ph) 

°/a  CHRX  Channelized  Receiver 

7a  out  =  chrx(inp,fs,Wc,fl,fh,Nz,dc,ph) 

7a 

7a  INPUTS 

7a  inp  -  input  signal 

7a  fs  sample  frequency  of  input  signal 

7a  Wc  -  channel  bandwidth 

7a  fl  -  lower  ChRx  frequency 

7a  fh  -  upper  ChRx  frequency 

7a  Nz  -  one-sided  zero  padding  length  (for  'filtfilt') 

7a  dc  -  downconvert?  (l==yes,0==no) 

7a  ph  -  downconversion  starting  phase 

1 

7.  OUTPUT 

7a  out  -  channelized  receiver  output  matrix 

7. 


7. 

7a  Brett  D.  Gronholz 

7a  EENG  799  —  Summer/Fall  2004 

7. 


if  nargin  ~=  8, 

error ('Not  enough  input  arguments!') 

end 

if  fs  <=  0, 

error( 'Sample  frequency  (fs)  must  be  positive  and  non-zero.') 

end 

if  Wc  <=  0, 

error ( 'Channel  bandwidth  (Wc)  must  be  positive  and  non-zero.') 

end 

if  fl  <=  0, 

error('Lower  frequency  (fl)  must  be  positive  and  non-zero.') 

end 

if  fh  <=  0, 

error('Upper  frequency  (fh)  must  be  positive  and  non-zero.') 

end 

if  mod( (fh-f 1) /Wc , 1)  ~=  0, 

error (' "fh-fl"  must  be  an  integer  multiple  of  "Wc"') 

end 

if  Nz  <  0, 
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error ('Zero  padding  length  (Nz)  must  be  positive  or  zero.') 

end 

if  (dc  ~=  1  &  dc  ~=  0) , 

error (' Invalid  parameter  (dc)') 

end 

N  =  4;  %  filter  order 

Nc  =  (fh-fl)/Wc;  number  of  channels 

delt  =  1/fs;  %  delta  t 

inp  =  [zeros (l,Nz)  inp  zeros (1 ,Nz) ] ;  /  zero-pad 

tdc  =  [0 : delt : (length(inp) *delt-delt)] ;  %  time  vector  for  downconversion 
for  i  =  l:Nc, 

if  dc  ==  1 ,  7.  downconvert 

[b,a]  =  butter (N,Wc/(f s/2) )  ;  7,  LPF  coeffs 

inpx  =  real (inp . *exp(-j*2*pi* ( (fl+(i-l) *Wc) *tdc+pi/180*ph) )) ; 
elseif  dc  ==  0,  7,  don't  downconvert 

[b,a]  =  butter (N,  [(fl+(i-l)*Wc)/(fs/2)  (fl+i*Wc)/(f s/2)] ) ;  7,  BPF  coeffs 

inpx  =  inp; 

end 

outt(i,:)  =  real  (f  iltf  ilt  (b ,  a,  inpx) )  ;  7  filter  input  sig 

out(i,:)  =  outt (i , Nz+1 : end-Nz)  ;  7  remove  zeros 

end 
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function  Cout  =  csm(x,Nfft,dim,clr) 
°/a  CSM  Cross-Spectral  Matrix 
7a  Cout  =  csm(x,Nfft,dim,clr) 


7a 

7a  INPUTS 
7a  X 

7a  Nfft 

7a  dim 

7a  clr 

7a 

7a  OUTPUTS 
7a  Cout 

7a 


input  matrix 

#  of  points  in  FFT  (to  form  SSM) 

dimension  (l==col-by- col ,  2==row-by- row) 

set  autocorrelation  terms  to  zero  if  1,  normal  if  0 


Cross-Spectral  Matrix  (CSM) 


1 

7a  Brett  D.  Gronholz 

7a  EENG  799  —  Summer/Fall  2004 

1 

if  nargin  ~=  4, 

error ('Not  enough  input  arguments!') 

end 

if  Nfft  <  1; 

error('Nfft  must  be  a  positive,  non-zero  integer.') 

end 

if  (dim  ~=  1  &  dim  ~=  0) , 

error (' Invalid  parameter  (dim)') 

end 

if  (clr  ~=  1  &  clr  ~=  0) , 

error (' Invalid  parameter  (clr)') 

end 

X  =  fft(x,Nfft  ,2)  ;  7a  FFT  of  rows  of  input  matrix 

ndim  =  size(X,mod(dim,2)+l) ;%  input  matrix  length  along  specified  dimension 
noth  =  size(X,dim);  7a  input  matrix  length  along  other  dimension 

clear  Cout 


if  dim  ==  1 , 

7.  Column-by-Column  Correlation 

Cout  =  zeros  (ndim,  ndim) ;  7a  initialize  matrix 
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%  C-by-C  correlation 


Cout  =  X’*X/noth; 
elseif  dim  ==  2, 

y„  Row-by-Row  Correlation 
Cout  =  zeros (ndim, ndim) ;  %  initialize  matrix 

Cout  =  X*X’/noth;  %  R-by-R  correlation 

end 

if  clr  ==  1 , 

°/„  Set  Diagonal  (Autocorrelation)  Elements  to  Zero 
for  i  =  l:ndim, 

Cout(i,i)  =  0; 

end 

end 
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function  Cout  =  ctm(x, Nfft, dim, clr) 
°/a  CTM  Cross-Temporal  Matrix 
7a  Cout  =  ctm(x,Nfft,dim,clr) 


7a 

7a  INPUTS 
7a  X 

7a  Nfft 

7a  dim 

7a  clr 

7a 

7a  OUTPUTS 
7a  Cout 

7a 


input  matrix 

#  of  points  in  IFFT  (to  form  TTM) 

dimension  (l==col-by- col ,  2==row-by- row) 

set  autocorrelation  terms  to  zero  if  1,  normal  if  0 


Cross-Temporal  Matrix  (CTM) 


1 

7a  Brett  D.  Gronholz 

7a  EENG  799  —  Summer/Fall  2004 

1 

if  nargin  ~=  4, 

error ('Not  enough  input  arguments!') 

end 

if  Nfft  <  1; 

error('Nfft  must  be  a  positive,  non-zero  integer.') 

end 

if  (dim  ~=  1  &  dim  ~=  0) , 

error (' Invalid  parameter  (dim)') 

end 

if  (clr  ~=  1  &  clr  ~=  0) , 

error (' Invalid  parameter  (clr)') 

end 

X  =  if  ft  (x,Nf  ft ,  1)  ;  7a  IFFT  of  columns  of  input  matrix 

ndim  =  size(X,mod(dim,2)+l) ;%  input  matrix  length  along  specified  dimension 
noth  =  size(X,dim);  7a  input  matrix  length  along  other  dimension 

clear  Cout 


if  dim  ==  1 , 

7.  Column-by-Column  Correlation 

Cout  =  zeros  (ndim,  ndim) ;  7a  initialize  matrix 
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%  C-by-C  correlation 


Cout  =  X’*X/noth; 
elseif  dim  ==  2, 

y„  Row-by-Row  Correlation 
Cout  =  zeros (ndim, ndim) ;  %  initialize  matrix 

Cout  =  X*X’/noth;  %  R-by-R  correlation 

end 

if  clr  ==  1 , 

°/„  Set  Diagonal  (Autocorrelation)  Elements  to  Zero 
for  i  =  l:ndim, 

Cout(i,i)  =  0; 

end 

end 
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